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Abstract 



We investigate the spatial statistics of the energy eigenfunctions on large quantum 
graphs. It has previously been conjectured that these should be described by a Gaus- 
sian Random Wave Model, by analogy with quantum chaotic systems, for which such 
a model was proposed by Berry in 1977. The autocorrelation functions we calculate 
^\ • for an individual quantum graph exhibit a universal component, which completely de- 

r) , termines a Gaussian Random Wave Model, and a system-dependent deviation. This 

deviation depends on the graph only through its underlying classical dynamics. Clas- 
sical criteria for quantum universality to be met asymptotically in the large graph limit 
(i.e. for the non-universal deviation to vanish) are then extracted. We use an exact field 
theoretic expression in terms of a variant of a supersymmetric cr model. A saddle-point 
analysis of this expression leads to the estimates. In particular, intensity correlations 
are used to discuss the possible equidistribution of the energy eigenfunctions in the 
\j~\ • large graph limit. When equidistribution is asymptotically realized, our theory predicts 

ps) \ a rate of convergence that is a significant refinement of previous estimates. The univer- 

^D ' sal and system-dependent components of intensity correlation functions are recovered 

by means of an exact trace formula which we analyse in the diagonal approximation, 
l/^ ■ drawing in this way a parallel between the field theory and semiclassics. Our results 

^P \ provide the first instance where an asymptotic Gaussian Random Wave Model has been 

established microscopically for eigenfunctions in a system with no disorder. 



Keywords : Quantum ergodicity. Random Wave Model, Criteria for universality. Rate 
of universality. Trace formulae. Nonlinear supersymmetric a model. 

1 Introduction 

Gaussian Random Wave Models are commonly used to describe the statistical proper- 
ties of the energy eigenfunctions of chaotic quantum systems. The original idea was 
introduced in 1977 by Berry IIT3JI , who proposed that a random function t// with Gaus- 
sian distribution 

could, in the semiclassical limit, reproduce all the spatial autocorrelation functions 

I r '' '' 



of a chaotic eigenfunction i^„ of energy e„. Here, 5 is a small volume that shrinks in the 
semiclassical limit but does so slowly enough to contain an increasing number of os- 
cillations of i^„, and/? in ([TJ is 1 if time-reversal symmetry is conserved, in which case 
xj/ is chosen real, and 2 if this symmetry is broken, in which case ij/ is complex. From 
a semiclassical calculation of C{ri,r2), Berry deduced that the covariance c{r\,r2) in 
([TJ is the free quantum propagator from r^ to r\. 

This is one of the central conjectures in the field of quantum chaos. Essentially, 
it asserts that the local statistics of quantum chaotic eigenfunctions correspond, in the 
semiclassical limit, to those of random superpositions of plane waves, and so are uni- 
versal. 

Following Berry, the universal Gaussian Random Wave Model has been refined to 
incorporate systems-specific features. For example, in quantum billiards, it does not 
fulfill the necessary boundary conditions. In this case, Hortikar and Srednicki 1291 
suggested replacing the covariance with the semiclassical approximation lfT4l |27l to 
the propagator of the system. This Gaussian model satisfies the boundary conditions 
and has the property that the direct path contribution to the semiclassical formula cor- 
responds to Berry's conjecture. Further understandings and refinements of this system- 
dependent Gaussian Random Wave model are given in B2l 14311441 . for example. 

It is important to emphasize that to-date effort has mainly been directed towards 
deriving the consequences of the Random Wave Model and its refinements, assum- 
ing its validity. Numerical tests strongly support the predictive value of the Random 
Wave Model. However, in no system has its validity yet been established or derived 
microscopically. 

We tackle here the problem of the validity of the Gaussian Random Wave Model 
on quantum graphs (a variant of this model was introduced in |f26l). Quantum graphs 
are favorable systems to gain some insights on the mechanisms responsible for random 
waves models to hold because, depending on their topology and their boundary condi- 
tions, their behaviors range from chaotic 1341 l35l . where a random model is expected 
to hold, to intermediate |I8||7||3T|[10||30], where such models should fail. Without any 
prior assumption on the nature of the quantum graph, one can evaluate its autocorrela- 
tion functions 



where {jc,),eN< and {j'yljeN, are points on the graph. In fact, we focus on the autocorre- 
lations for q - p - \, and for q - p with {jc,-),eN = {jylyeN ■ The other autocorrelations 
are believed to vanish due to additional complex phases that fluctuate strongly. The 
result obtained for q — p - 1 is exact and yields a universal covariance c which de- 
fines the unique candidate for the Gaussian model on quantum graphs. It should be 
emphasized that this does not contradict the construction of Gaussian Random Wave 
Models with a system-dependent correction in analogy to Urbina and Richter's guess 
for billiards B2l . Indeed, our autocorrelation functions are defined by averaging over 
the whole energy spectrum. Such an average on Urbina and Richter's random functions 
also kills the system-dependent correction and leads to a covariance given by the free 
propagator, namely, to Berry's universal model. The system dependency found in the 



autocorrelation functions of higher degree evaluated here is of different nature. It is not 
a refinement of a universal Random Wave Model, but it rather measures how chaotic a 
given quantum graph is from the energy eigenfunctions perspective. Interestingly, this 
non-universal term is found to depend on the quantum graph only through its classical 
dynamics. This provides us a way to estimate the deviation from quantum universality 
in terms of a classical quantity, and so to discuss criteria for the Random Wave Model 
to hold in the large-graph limit, and cases where it fails, such as Neumann star graphs 
IfTTI . Our results provide the first instance where an asymptotic Gaussian Random 
Wave Model has been established microscopically for eigenfunctions in a system with 
no disorder 

It would be a major achievement to show that our result for the autocorrelation 
functions ((Sj in the case of quantum graphs also applies to other quantum systems. If 
this is the case, the deviations from universality vanish in the semiclassical limit in 
chaotic billiards, which explains why such deviations have indeed never been found, 
whereas they must prevail over the universal part in non-chaotic systems. For chaotic 
systems, the corrections would reveal the rate of approach to universality as h —> 0. 
Finally, if such a formula was found, its ability to describe systems with mixed phase 
spaces could be studied and compared with the empirical results ||2] |3| and alternative 
approaches based on bifurcation theory and singularity-dominated strong fluctuations 

ma. 

The moments and autocorrelations of second degree (i.e. intensity correlations) 
play a particularly important role in quantum chaos, because they suffice to measure 
the spreading of the energy eigenfunctions, and they can be rigorously controlled. Ac- 
cording to |40|, the high energy eigenfunctions of a classically ergodic system should 
become uniformly spread over the surface of constant energy, a property known as 
quantum ergodicity. This claim has found rigorous proofs in ifTTl . lfT6l and ll50l for 
example, where the authors consider compact manifolds with ergodic geodesic flows, 
quantized ergodic maps and ergodic billiards respectively. The main tool used in these 
works is an Egorov estimate, which, in the case of quantum maps, reads 

WuJ^OpifW'l, - Op(J o m')\\ < const ■ h, (4) 

for M a map, / any smooth function on the configuration space, and U and Op{f) their 
quantized analogs. A version of (|4]i also holds for continuous Hamiltonian systems. 
However, the Egorov method does not provide any information on the rate with which 
quantum ergodicity is reached. This much harder problem is investigated in 1211 |47l 
|48l|49l[Tl[T][33. 

In fact, quantum ergodicity is significantly more difficult to tackle on quantum 
graphs than on other chaotic systems. The reason is the non-existence of a determin- 
istic classical map, and hence, of an Egorov estimate. In |9l, quantum ergodicity is 
proved for graphs related to quantum maps by using the Egorov property on the under- 
lying quantum maps. On the other hand, it is shown in |[10l [11] [30l that some graphs, 
namely star graphs, are not quantum ergodic. Here, our result for the autocorrelation 
functions (O with p -2 enables us to expound a criterion for graphs to become quan- 
tum ergodic. A summary of our results in this special case has already been given 
in ll24l . Moreover, our method also yields the rate of quantum ergodicity in terms of 



the classical dynamics, when quantum ergodicity does occur. The result obtained is a 
significant refinement of the previous estimates in ifTSl . 

The reader not interested in the derivation of the formulae can directly jump to Sec- 
tion]?] where the final formulae are given and exploited. The rest of the text is structured 
as follows. In the sections [2] and [3] quantum graphs are defined, and the autocorrela- 
tion functions together with other statistical quantities of interest are introduced. In 
particular, a first type of trace formulae is developed in 13.41 and 13.51 An exact field 
theoretic expression for the autocorrelation functions is developed in Section [H and a 
second type of trace formulae is presented in 14.21 Then, two different contributions 
to the exact expression for the autocorrelation functions are extracted and calculated 
in sections |5] and |6l Section [7] compares these two contributions and illustrates them 
with a few examples. Section |8] discusses our results and gives an outlook on possible 
implications. 

2 Quantum Graphs 
2.1 Definitions 

A metric graph G is a set of V € N points, called the vertices, and of B e N bonds 
of positive lengths L = (Li , ■ ■ ■ , Lb) linking some pairs of vertices. The topology of a 
graph is determined by its connectivity matrix C, namely the V x V matrix 

Cij - Cjj = #{bonds connecting the vertices / and j}. (5) 

If Cij - and C/j < 1 for all i, j € Ny, the graph is said to be simple. The valency v,- 
of a vertex / 6 Ny is defined by v, = Tj]=\ Cij- The valencies are always all supposed 
positive. A point on a graph is specified by a pair {b, Xb), where b eNg determines the 
bond and xt, 6 [0, Lh] determines the position of this point on b. 

Each bond of a metric graph can be traversed in two possible directions, denoted 
by d e {+, -). A pair/? - (b, d) then denotes a directed bond, and p = (b, —d) stands 
for its reverse partner. The vertex from which a directed bond /? emerges is written ofi 
and the vertex to which it leads is written tp. In particular, o/3 - tp is always fulfilled. 
We suppose the set of directed bonds to be ordered so that, by abuse of language, any 
directed bond^S can also be seen as an element of Nib. 

A quantum graph is a metric graph G that is turned into a quantum system. In order 
to do this, the C-linear space 

B 

-7/ = {^P = ^„\^,„ ^l ^l 6 l\\Q, Lb\)] (6) 

b=\ 

is introduced, and its elements are referred to as wave functions. This space is endowed 
with the scalar product defined by 



cF,(i)) = y, f V;w^ft(^¥^ (7) 



f^-'o 



for any *?, O € 1-1. The number if>b(xb) is interpreted as the value of the wave function 
*!* at the point (b, x/,) of G. One can define an operator H acting on '7/ as 

B B 

//0^/,-0-'a;- (8) 

b=l b=l 

This is the expression of the free quantum particle Hamiltonian on each bond. The 
restriction of H on the subset "Ko c '7Y of wave functions vanishing at the vertices is 
symmetric. A wave function ^ € 'Hq is called a Dirichlet wave function. A Schrodinger 
operator on a metric graph (and thus a quantum graph) can be defined as a self-adjoint 
extension of H. However, we will follow a slightly different definition using the scat- 
tering approach 1341 . We first give a brief overview of this approach and then discuss 
its relation to self-adjoint extensions of H. 

For any real number ^ > 0, the solutions of the equation //*!' - k^^' form the 
subspace 



fl/j€C, V/3 = (^,£/)€N2b[, (9) 



[ b=l d=+- 

where, for /? e Ng and d e {+, -), 

-Cbdik) = e^'i''-"-^). (10) 

A wave function in ^{k) is then characterized by 2B waves of wave number k, each 
of which carries a complex amplitude ap corresponding to its value at the mid-point of 
the bond. 

Let us introduce 2B formal symbols le^j), fi € N2B, and the set J{. of their possi- 
ble linear combinations over C The set y[ is a 2B-dimensional C-linear space called 
amplitude space, and it is endowed with the hermitian scalar product defined by 

<e/)'k/s> = (5'^,/;'- (11) 

It can be seen as the direct product J[ - J[b ® J[d of a B-dimensional bond space Jib 
and a two-dimensional direction space yi,;. For each A: > 0, there is a natural one-to-one 
mapping 

B IB 

^^®Yj ""'d'^hdik) ^ la) = ^ aplep) (12) 

h=l d=+- /3=1 

between ^(A) andyi. If ^1 h-» \a{) and ^2 •-» l«2) by this mapping, the scalar products 
in the spaces ^ik) and J{ translate 



(>I'i,»I'2)= ai 



ai), (13) 



where (f[ stands for the first Pauli matrix acting on Jij, and L denotes the 2B x 2B 
diagonal matrix 

V/3 = ^isfi'h- (14) 



Here and henceforth, the length of a directed bond is the length of the bond on which 
it is supported. In particular, Lp - La is always fulfilled. The identity ( fT3] l shows that 
the mapping (fTZb does not preserve length and orthogonality in general. 

In the scattering approach to quantum graphs the values at each vertex / G Ny of 
the V, waves emerging from this vertex and of the v, waves incoming to this vertex 
are related through some fixed matrix cr' . If Iflom) and \a'- ) denote the v, -dimensional 
vectors containing the values at vertex / of the emerging waves and of the incoming 
waves respectively, this relation reads 

l«U> - ^'lO- (15) 

A wave function ^P e ^{k) conserves the probability current if and only if the V matri- 
ces cr' are all unitary. The components of the V outgoing and incoming vectors \a'^^^ 
and \a' ) can then be grouped together to form the 2B-dimensional vectors lOout) and 
lain) respectively. These vectors are related to \a) in (fT2b through 

Kut) = r'-«|a> and \a;,) ^ T (k)\a) (16) 

where T{k) is the 2B x 2B diagonal matrix T{k) — e''^'- . This matrix contains the phases 
gained by the 2B waves of wave number k when they travel along half the bonds on 
which they are supported. It is referred to as the propagation matrix. Moreover, the V 
identities (flST l become 

Iflout) = 5|ai„), (17) 

where S is the 2B x 2B unitary matrix, called scattering matrix, defined by 

V.-l?''^ ''f-'P-' (18) 

'^'^ \ otherwise 

Putting ( fTSI l and ( fTTl i together yields 

U{k)\a)^\a), with Uik) ^T{k)ST(k). (19) 

The 2B X 2B matrix U(k) is called the quantum map or evolution map of the graph. It 
is unitary since both T{k) and S are unitary. 

Equation ( fT9] l shows that imposing the conservation of probability current through 
fixed unitary matrices cr' restricts the possible amplitudes \a) and the possible wave 
numbers A: > 0. Indeed, the secular equation 

det(l-[/(A:)) = (20) 

must be satisfied for ( fT9b to admit non-trivial solutions. This equation is satisfied for a 
sequence 

< ^1 <k2 < ■ ■■ <ky < ky+i < . . . ^ oo (21) 

called the spectrum of the quantum graph, and the square of these wave numbers are 
the quantized energies. If the bond lengths L\,. . .,Lb are independent over Q, there is 
typically a normalized vector \a^) in J{ for any v € N that satisfies U{kyy\d^) - |a'') and 
so that any other vector satisfying this equation is of the form z\d') for some z € C The 



vector la'') then provides the amplitudes of the eigenfunction ^f satisfying //*!"' = k^^''' 
by the mapping ( fT2] i. IncommensurabiUty of the bond lengths and this non-degeneracy 
property will be assumed henceforth. 

It is well-known IZSl that the mean number of allowed wave numbers in [0, K] is 
A^(^) = Kd, where the mean level density d reads 

^.|. (22) 

For any k > 0, the unitarity of U(k) ensures the existence of an orthonormal basis 
{\n, ^))«eN28 of C^* and of 2B real numbers {(^„(^))nENo« such that 

f/(/t)|«,;t) = e'*W|n,yt). (23) 

These sets can be ordered by imposing the inequalities 

-2n< 02b(O) < (piB-iiO) <...< 02(0) < <^i(0) < (24) 

and by requiring the 2B eigencurves k i-> (p„{k) to be C°°. This smoothness condition 
can indeed be realized since the map U{k) depends on k in an analytic way. Taking a 
derivative with respect to k on both sides of ( l23b leads to 

(f)'„{k) - <«, k\L\n, k) e [Lmin, imax], (25) 

where £„,!„ and L^ax denote the minimal and maximal bond lengths on the graph. 

A quantum graph is time-reversal invariant if its quantum map satisfies tr ( U{k)'^) - 
tr U(k)" for all A; > and integers n. Here and henceforth, the generaUzed transposition 
A^ of a linear transformation A is defined by 

A^ = crfA^trf, (26) 

A^ being the transpose of A. It satisfies A'^ — A. Since Tik)"^ - T(k), a graph is 
time-reversal invariant if and only if its scattering matrix satisfies tr (5''^] = tr 5" for 
all integers n. Obviously, S''' = S implies time-reversal invariance. Note, however that 
replacing 

S ^S' ^ e-'^Se'^ \a) ^ \a'} = e^'V) (27) 

where 9 = diag(0i, . . . , 02b) is a diagonal real matrix is equivalent to choosing a differ- 
ent reference phase for the amplitudes. We will call such a transformation a (passive) 
gauge transformation - it neither affects the spectrum nor the condition described above 
for time-reversal invariance. The latter can now be reformulated: a quantum graph is 
time-reversal invariant if and only if there is a (possibly trivial) gauge transformation 
S 1-^ S' - e^'^Se'^ such that 5"^^ = S'. For time-reversal invariant graphs we will 
henceforth assume that the reference phases have been chosen such that 5^ = 5' holds. 
There remains a residual gauge freedom to which we will return later when we discuss 
the wave function statistics in quantum graphs. 

The set of all time-reversal invariant graphs form the orthogonal symmetry class, 
and the set of all quantum graphs violating this property form the unitary symmetry 



class. We will frequently use the parameter k which takes the values 

1 in the unitary class, and 



(28) 
2 in the orthogonal class. 

Note that the parameter k that we use here is linked to the parameter yS used in random- 
matrix theory to distinguish symmetry classes by a- = Ijfi. 

We have already mentioned that the scattering approach described above is not the 
only way to define quantum graphs. The other frequently used definition is based on 
self-adjoint extensions of // in (O defined on the Dirichlet domain "Ho (see jTSl and 
references therein). A complete description of all possible self-adjoint extensions was 
given in Il33ll . In general each self-adjoint extension is equivalent to energy-dependent 
matrices cr^^''{k) relating the outgoing amplitudes to the incoming amplitudes of 'I' e 
^{k) at each vertex / instead of ( fTSl ). These matrices can then be grouped together to 
form a global unitary scattering matrix S^^{k) as in ( fTSl l. and a global quantum map 
U^^(k) = T{k)S^^(k)T{k) satisfying the secular equation ( |20| ). The two definitions of 
quantum graphs have a certain overlap as there is a subset of self-adjoint extensions 
which leads to energy-independent scattering matrices. It is shown in |6l and ifTSl 
that any scattering matrix S ^^ defining a self-adjoint operator H admits a limit S ^^ 
as k tends to infinity, and moreover, it is argued in ||6l that a scattering matrix S^^ 
and its limit S^^ share the same spectral statistics. The coincidence of these statistics 
comes from the fact that they are properties at asymptotically large wave number k. 
Hence, one can deduce that the eigenfunction statistics of 5*^^ and 5^^ also coincide. 
As a consequence the eigenfunction statistics of quantum graphs defined following the 
self-adjoint extension approach can be recovered from the eigenfunction statistics of 
quantum graphs defined through the scattering approach by substituting S^^ for 5''^^. 

Henceforth, the scattering matrix S always refers to the matrix in ( fTsT i obtained 
from the scattering approach. It can be any 2B x 2B unitary matrix such that S^'/j 
vanishes if t/3 + op . A possible choice is the so-called Neumann scattering matrix, 
which is defined at each vertex ; e Ny by 

4/^=--^A/3'' W^o/i'^i. (29) 

Quantum graphs with this choice of scattering matrix at each vertex will be called 
Neumann quantum graphs. 

In general, a quantum graph is then specified by a pair (G, S ) where G is a metric 
graph and 5 is a scattering matrix on G. The class of possible scattering matrices S 
on G contains all the asymptotic matrices S ^^ obtained from the self-adjoint extension 
approach. There are however some scattering matrices that are acceptable from the 
scattering point of view but not from the second approach. An example is given by 
the Direct Fourier Transform (DFT) graphs f25l, for which the scattering processes at 
vertex / e Ny are described by the v,- x v; unitary matrix 



cr^,^ = -^e-""^^i~, Mtp = op' = U (30) 



where n' is a surjective assignment of an integer in N,., to each directed bond around 
i such that n'ifi) - n'(J3). With these matching conditions, the wave functions {*I"')veN 



obtained from the amplitudes {|a''))v£N and the spectrum {fcvJveN by ( fT2] l are in general 
not orthogonal to each other in '7/, which shows that H acting on the wave functions 
satisfying ( [30l l is not self-adjoint. By contrast, the Neumann scattering matrices ( [29] l 
do lead to a self-adjoint Laplace operator 

Both examples of scattering matrices were defined in terms of symmetric unitary 
matrices at each vertex cr' - cr' . As a consequence S"^ - S and the a quantum graph 
obeys time-reversal symmetry. One may break time-reversal symmetry by adding a 
magnetic field to the graph. In the scattering approach adding a magnetic field which 
is constant on every bond is straightforward. Let A be the diagonal matrix that contains 
the magnetic field strengths. It obeys Ap - -An. The corresponding quantum map is 

U{k) = e'(*+'4W25e'(*+'4W2 = T{k)SAT(k) (31) 

and the magnetic field effectivelyjust changes the scattering matrix 5 >-> Sa - e'^^^^Se'^^^^. 
If 5 - S"^ and the graph is multiply connected (that is, it contains cycles) then the mag- 
netic field generally breaks the time reversal invariance. 

Henceforth, the metric graphs G considered are assumed simple. The reason for 
this assumption is to simplify some notations and calculations. However, if a graph 
contains a directed bond fi such that of} - tp, a Neumann vertex can be added on the 
bond b supporting /3 to destroy the loop b without modifying the quantum dynamics. 
Similarly, if the graph has two directed bonds /3,/3' such that o/3 - o/3' and t/3 - t/3', 
a Neumann vertex can be added on the bond b supporting /3 to destroy this parallel 
connection without modifying the dynamics. Hence, any graph can be made simple 
by adding sufficiently many Neumann vertices, and this process does not change the 
quantum dynamics. One can thus assume the graph simple without loss of generality. 

2.2 Classical Dynamics 

With any quantum graph, one can associate a bistochastic classical map M defined by 

M^^< = Wi^flikt = 15^^-1', (32) 

where U{k) is the quantum map and S is the scattering matrix. The matrix M describes 
a Markov process on the graph, which is the classical counterpart of the quantum dy- 
namics defined by S . The uniform vector 

1 '^ 
11). ^gM (33, 

is an eigenvector of M of eigenvalue 1, and its hermitian conjugate (1| is a left eigen- 
vectors of M of eigenvalue 1 . Besides, the Perron-Frobenius theorem Ii28l ensures that 
the spectrum of M Ues on or within the complex unit disc. 

A graph is said to be ergodic if and only if, for any /3,/3' € N2B, there is a discrete 
time « 6 N for which the transition probability {epi\M"\ep) is positive. This condition 
is equivalent to the non-degeneracy of the eigenvalue 1 of M. Any non-ergodic graph 
(G,S) is the union of several ergodic components, that is (G,S) - Ut=i(^"'^') f°i' 



some integer A: > 1 . The eigenvalue 1 has degeneracy k, and the k vectors that are uni- 
form on one component (G,, 5,) and zero on the others form a basis of this eigenspace. 
Let us write M^ = e^^'^M for an ergodic classical map M and for some e > 0. The 
sum of all classical paths fromyS e Nib to /?' e N2B followed with M^ can be written 

It becomes singular as e approaches zero due to the eigenvalue 1 of M. Let M - 
Dm + Nm be the Jordan decomposition of M into a diagonalizable part Dm and a 
nilpotent part Nm commuting with each other Let {/lj)jeN,i, be the 2B eigenvalues of 
Dm, and let {|y))jeN,B be corresponding normalized eigenvectors in J[ with |1) as in 
( |33] |. Then, it is straight forward to check that (l|j) - 61 j. This fact enables one to 
extract the singular part of ( l34l i and write 

^ = ^— ^IIXIK/^. (35) 

where the remainder R^ is such that R = lim^^o/^f exists and satisfies {\\R - and 
7?|1) = 0. The first and second terms in the right-hand side of ( [35] l will respectively be 
referred to as uniform and massive components. If m, = 1 - /I,- for i -2, . . . ,2B denote 
the 2B - I non-zero eigenvalues of 1 - M, the massive component satisfies 



XrR = Y, 



2B , 

'-^. (36) 



The eigenvalues {m,},EN,B of 1 - M are called masses. They all lie in the closed disc 
of radius 1 and centered at 1 in the complex plane, and the zero mass nii = is non- 
degenerate. 

3 Eigenfunction Statistics 

3.1 Wave function correlation functions 

Let {G,S) be a quantum graph, [ky] be its spectrum, and {a") c C^* be a set of nor- 
malized amplitude vectors defining the eigenfunctions {*!"') as in (fT2] i. Let us consider 
2B complex random variables ap and investigate the existence of a joint probability 
density function (f{a) - ip{a\, . . ., 023) satisfying 



> A-=0 1=0 I ^ ' k,<K ^ '^ k=Q 1=0 

= I f]4f]fl^;^(«)flfaVa, (38) 



k=0 1=0 



for any choice of ySo, . . . ,y6,,-i ,/?,',, . . . ,j6' 1 6 ^2b with p,q 6 Nq. Here, the measure 
da* da denotes the product of the 2B flat Lebesgue measures dalda/j in the complex 
plane, and the notation {0}v for a2B x2B matrix O stands for {0}v = (a^lOla"). 
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The first line in (|37] | defines the wave function correlation functions. The peculiar 
factor T^TTT- in this definition is introduced for further calculational convenience. For 
large graphs with extended wave functions this factor is expected to be close to unity. 
Indeed, it has generally a tiny effect on the wave function statistics. It will be seen later 
that with the inclusion of this factor ( [37] i does not depend on the particular values of 
the incommensurate bond lengths. Moreover, performing an average over the spectrum 
of the quantum graph in presence of this factor, such as in ( |37] |. amounts to averaging 
the same quantity over all the eigenfunctions |«, k) of U{k) and then integrating over all 
k G (0, oo). Indeed, it is proven in ITZl that graphs with incommensurate bond lengths 
obey 

1 trL 1 r^ 1 "^ 

lim y {0)t = lim - —y{n,k\0\n,k}'^dk. (39) 



N{K) ^^ 2B{L\. 



n=\ 



for any 2B x 2B matrix O and any non-negative integer q. 

The identity (|39] l shows that the joint probability density function ip{a) in ( l38l l is 
normaUzed. Indeed, choosing ^ = in this formula leads to 



(D.nm^y^^^limir^Vl.l. 
K^^ N(K) f-i 2B{L}y K^^ K Jo 2B ^ 



IB 

(40) 

ky<K ' 



Moreover, it also provides an exact expression for the covariance of ip{d). Indeed, the 
equality ( |39] | with q-\ and O - \ep'){ep\ yields 



1 <-^r 

(al,ag) = lim > (>'k«')(efl|v) 



^, (41) 



K^c. N{K) ^^ 2B{L) 

1 r'^ 1 "* 

= Um — I — > {n,k\eR'){ea\n,k)dk 

K^oo K Jo 2B j^^ ^ '^ '^ '^^ ' 2B 

by orthonormality of the families {\n, k)]nen^B- This derivation of the covariance relies 
on the incommensurability of bond lengths. We will show later in Subsection l3.4l that 
the restriction to incommensurable bond lengths can be lifted. 

Further properties of the joint probability density can be derived considering its 
invariance under gauge transformation of the form described in ( |27] i. This discussion 
has to treat systems with and without time-reversal invariance separately and we will 
start with the unitary class (broken time-reversal invariance). In this class we are free 
to choose a gauge and one expects that all correlation functions which are not gauge in- 
variant will vanish. This implies that the non-trivial correlation functions (|37] | have the 
same number of complex conjugated amplitudes as non-conjugated amplitudes (that is 
p — q). It is thus suflicient to consider the autocorrelation functions 

C[„] = (k|^..K,_,|2^. (42) 

where [or] = [ao, . . . , ctq-i] is a vector containing ^ e No directed bonds aj e N2b- The 

integer q is called degree of C[„] . 

Of particular interest to us are the moments 

Ma,c, = C[,xa] = (iflj'") (43) 
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and the first non-trivial autocorrelation functions 

Caa' = C[aa'] = {\aa\^\aa'\^) (44) 

which form the symmetric intensity correlation matrix Caa' = Ca'a- 

For time-reversal invariant systems the property U(k)'^ - U(k) of the quantum map 
implies that one may always choose the phase of its eigenvectors \n, k) - J/j fl„,^(A;)|e/3) 
such that ai,fi{k) - a^d,k)* . An equivalent statement is that the wave function on the 
graph can be chosen real. This has strong implications on wave function statistics - for 
instance the autocorrelation functions C[„] defined in ( l42l i are invariant under replacing 
any directed bond in cc = {oq . . . ,aq-\) by its reverse partner a, i-> a,. The joint 
probability density function then reduces to a product 

ip{a) = 5*(a+ - a*_) ipKi(a+) (45) 

where a+ (a) is the B-dimensional vector containing the amplitudes for directed bond 
a - (b, d) with positive (negative) direction index d. For a quantum graph in the 
orthogonal class it is thus sufficient to consider only the correlation functions in ([37j 
for which all directed bond have a positive direction index. As in the unitary case we 
also expect for the orthogonal case that correlation functions that depend on a local 
gauge vanish exactly. Note, that in the orthogonal case not all gauge transformations 
dZTJ i are allowed. In order to preserve the properties 5^ = 5 and On^ik) - a^g{ky 
only gauge transformations with fl/j - -9a are allowed. Again the only non-trivial 
correlation functions are the autocorrelation functions (l42l i. 

3.2 Circular and Gaussian Random Waves Models 

For a large well-connected quantum graph in the unitary symmetry class the quantum 
map U{k) does generally not have any symmetries. Moreover, in a complex network 
(e.g. a randomly chosen connected graph) the neighborhood of any bond looks statis- 
tically the same. By analogy with the circular ensembles of random matrix theory one 
is inclined to guess that the joint probability density function (p(a) for the eigenvectors 
of the 2B X 2B matrix U{k) defined in ( l38l l is invariant under transformations a >-^ ua 
for unitary matrices u. This implies that the vectors a are uniformly distributed over 
the unit sphere in C^^. We will call the guess 

^cv(a)=^^^^6(l-\\af) (46) 

the Circular Random Wave Model for quantum graphs in the unitary class. The mo- 
ments and the intensity correlation matrix predicted by the Circular Random Wave 
Model read 



„„,„, =,ga^L^ =^(,_%^,««-.) 



(2B-H^-1)! (2B)?\ 4B 



^^■''•"' -(2B)(2B+1) -^W\~^^^^^'^ 



Yl 



(47) 







Figure 1 : Numerically evaluated intensity correlation matrices Caa' defined in (l44l) for 
a complete graph with V = 16 vertices and B = 120 bonds. In the upper two panels 
the DFT scattering matrices have been used at each vertex. The lower two panels are 
for Neumann scattering matrices. For the right two panels time-reversal symmetry 
has been broken by adding a magnetic field. The directed bonds a = (b, d) have been 
ordered as ((1, -),(2, -), . . .,(B, -), (1,+), . . .,{B, +)). For the graphs in the orthogonal 
class on the left side there are four identical blocks as Caa' - C&^a' - Caa' - C&^> . In 
the unitary case, note that the correlation matrix on the off-diagonal a = a' remains 
strongly peaked for Neumann scattering matrices. However the four blocks are no 
longer the identical (this is not obvious from the picture). For DFT scattering matrices 
the strong off-diagonal peak almost disappears in the presence of a magnetic field. 
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In the limit B ^ co (and constant degree q) one may replace the Circular Random 
Wave Model by the Gaussian Random Wave Model with the joint probability density 

^Gu(«) ^ ^ e-^*"^ (48) 

The predictions for the moments and the intensity correlation matrix in the Gaussian 
Random Wave Model read 

q\ 1 + 6aa' 

MGU,a,g = J^ and CGV,a.a' = ^^gy_ (49) 

which is equivalent to the leading order of the predictions ( |47l ) of the Circular Random 
Wave Model as B — > oo. 

If the eigenfunction statistics dJTl i of a family of quantum graphs in the unitary 
symmetry class are well reproduced by (pc\j(a) in ( |46] | or by (pcuia) in (l48l l. these 
formulae provide us with a universal Circular or Gaussian Random Wave Model, which 
gives access to all the statistical properties of the eigenfunctions. Notice that the exact 
calculation (HTt asserts that ( |48] | is the only possible Gaussian joint probability density 
function of the type ( l38l l. and hence, a non-universal Gaussian model cannot be realized 
on quantum graphs. 

Establishing the possible validity of the Gaussian Random Wave Model (|48] | would 
require the calculation of ( |37] i for arbitrary products of amplitudes a^. Note that the 
Gaussian Random Wave Model is consistent with the gauge principle, i.e. its prediction 
for any correlation function that is not explicitly gauge invariant vanishes identically. 
In what follows we will mainly focus on the explicitly gauge invariant autocorrelation 
functions (|42] |. However in subsection 14.21 we will show that some low order correla- 
tion functions that are not explicitly gauge invariant indeed vanish on the level of the 
diagonal approximation. 

When time-reversal symmetry is conserved one has to take into account that the 
amplitudes of counter propagating waves on the same bond are complex conjugates, so 
that the wave function is real. We can thus only expect that a universal joint probability 
function is invariant under a i-> ua where m is a unitary 2B x 28 matrix that respects 
reality of the wave function or, equivalently, that S h^ u^Su conserves S''~ - S . Such 
unitary matrices have the block structure 

r ""-) (50) 

in terms of the direction index d. Here u++ and m+_ are two B x B matrices which 
are only constrained by unitarity of u. Unitary matrices with this block structure obey 
m"^ - m' = ir' and are thus in fact orthogonal matrices with respect to 7"-transposition. 
The Circular Random Wave Model for the orthogonal class 

2^->(B-l)!,B,„ „,_/! 



^co(a) = ^^^^5«(a+ - a*) 6l- - ||a+|H (51) 

is the unique model which respects a+ = a* , the normalization \\a\\^ = 1 = 2||a+|p, 
and is invariant under the generalized orthogonal transformations ( fSOl l. It gives the 
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predictions 



^!(B-1)! _ q\ (^ q(q-\) , ^^, ^_^^ 

(52) 



Mco.a., - s: /':., = 7:^(1 -^^^^TTe+c?(B-') 



2i(B + q-iy. (28)1 \ IB 
Cco,«,.' = 4fi(fi+i) = ^^)2 [^--^+0{B ) 

The only difference in the leading order for large graphs is the term 6adi which ensures 
that the intensity correlation matrix is invariant under a i-> a. Note that the deviations 
in the next order are twice as large in the orthogonal case. 

In the limit B — > 00 one may again replace the Circular Random Wave Model by a 
Gaussian Random Wave Model with the joint probabiUty density 

^Go(«)^^A«.-«*).-'*^"' (53) 

where only one half of the coefficients is taken from a Gaussian ensemble while the 
other half remains fixed by the symmetry constraints. The moments and the intensity 
correlation matrix in this Gaussian Random Wave Model are just the leading order 
terms from ( 



Mco,y.q = ^2^ and Cco^a.a' = -^2 • (54) 

Note that the unitary and orthogonal universal Gaussian Random Wave Models 
( |48] l and (|53T l do not obey the normalization condition ||a|p - 1. In fact one has 

(ii«ii')gu = (ii«iOgo-i- (55) 

only as an average property while the variances 

are positive. Similarly, |flj^| cannot exceed one while the Gaussian Random Wave Mod- 
els have a finite probability for this event. The Circular Random Wave Models take all 
these constraints into account correctly. 

There is another obstruction to all the Random Waves Models (l46b . (l48T l. (ISTl i and 
dSSl l. The matching conditions at vertex / impose some correlation between the ampli- 
tudes supported on the neighboring bonds. This type of local and system-dependent 
correlations is ignored in the universal Random Wave Models. The most striking ex- 
ample consists in adding a Neumann vertex on some bond b of an ergodic graph. By 
doing so, the bond b is split into two new bonds b\ and Z72, which can be oriented such 
that (^1,-1-) — > (b2,+)- Then, the Neumann condition imposes |fli,+|^ - |fli,+P and 
l^fci-l^ = |flfeT-P- These strong correlations contradict the predictions WT\ and ( |52| |. 
Hence, a necessary condition for the universal Gaussian models ( l48l l and ( |53] | to be 
fulfilled in the limit of large graphs is that all the valencies tend to infinity. 

For a finite graph one should expect that none of these models reproduces the exact 
correlation functions. Indeed, any numerical evaluation of the wave function statistics 



15 



shows (amongst other things) an intensity fluctuation matrix that is far less uniform 
than the predictions from the Random Wave Models (see Figure[Tli. The deviations can 
only be expected to vanish as B — > oo and if certain other conditions that we are going 
to derive are also satisfied. 

3.3 Asymptotic Quantum Ergodicity 

Let G be a metric graph with B bonds. An observable on G is a family 

V^[VbeC\[0,Lb])\beNB} (57) 

of B real functions Vb{x) defined on the bonds of G. The mean value V of an observable 
V is defined by 

Notice that ^ - jl lis the volume of G. If an observable V is constant on each bond, 
one can simply write V - {Vb)be¥iB ^i^h Vb e K. The mean value of such an observable 
reads 

and is invariant under a global scaling of the bond lengths. 

Suppose now that S e U{2B) is a scattering matrix on G. The quantum graph 
(G, 5 ) is said to be quantum ergodic if and only if there exists a subsequence / i-» v(i) 

of density 1 such that 

(\iiv(i) yxuv(')) 
lim ^-—^, — ^ = V (60) 



?— >oo 



for any observable V. In this definition, ^'^ = ®b=\ ^l denotes an eigenfunction of H 
of eigenvalue k^. By assumption, it is unique up to multiplication by complex numbers 
(or by real numbers for the orthogonal class). 

The left-hand side of ( l60l l represents the mean value of the observable V in the 
eigenstate ^'^^'\ A straightforward calculation shows that 

-Li 



CV^ V^'') = y (\alf + \al_f) f ' Vb(x)dx 

ti Jo 

+2% y al*_al^ f " n(x)e2''^"(^-^)rfx (61) 

for the wave function *P'' with wave number ^v > and amplitudes fl]J^ and a^ as in 
(fT2l i. Since the observable V is assumed continuous on each bond, and since |a]^'^aJJ^| < 
1, the second term in the right-hand side of ( l6Tl l is 0{k^^). In the high energy limit this 
second term gives no contribution to the left-hand side of (|60] |. Moreover, the first term 
in the right-hand side of dMT l remains unchanged if the observable V is replaced with 
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the observable W defined by Wb = L^^ L'' Vbix)dx. These two remarks imply that, in 
the definition ( |60] | of quantum ergodicity, it is sufficient to consider observables that 
are constant on each bond, and this will always be the case in what follows. 

If the equality ( l60l l holds for any observable of vanishing mean V = 0, then it also 
holds for any observable W. In order to see this, it is sufficient to observe that W - tV 
has vanishing mean and to apply (|60] | to this new observable. Hence, without loss of 
generality, one can also restrict attention to observables V with V = 0. 

If the identity (|60] | is satisfied for any subsequence of eigenfunctions, the quantum 
graph is said to be quantum unique ergodic. In 1381 . it is shown that many short closed 
cycles, like the triangle/?! -^ Pi —^ f^s —^ /3i for instance, support eigenfunctions with 
arbitrarily high energies. These eigenfunctions, called scars, break quantum unique 
ergodicity. While these scarred eigenfunctions were obtained explicitly for Neumann 
quantum graphs, quantum unique ergodicity should certainly not be expected to hold 
on general finite quantum graphs. 

Moreover, quantum ergodicity is generally not realized on a finite quantum graph 
as well. This notion has thus to be replaced with a weaker one which we call asymp- 
totic quantum ergodicity. Let us consider an infinite sequence {(G/, S /))/eN of quantum 
graphs with increasing number of bonds Bi < Bi+\. We also suppose that the bonds of 
any G/ have bond lengths that satisfy 

Li 6 [Linin,imax] where < Lmin < imax < '>« (62) 

are independent of I. Such a sequence will be called increasing. We always assume that 
either all the graphs (G/, S i) are time-reversal invariant, or they all break this symmetry. 
The eigenfunctions of (G;, S /) are denoted by ^[, and similarly, all the quantities intro- 
duced above are indexed by I. Besides, a sequence {V/j/eN, where Vi is an observable 
on G/, is said to be acceptable if and only if the two conditions 

lim/^co V; = Voo exists, 

o<|y,,ii<yn,ax ^ ' 

are fulfilled. Then, an increasing sequence {{Gi,S i)]ien of quantum graphs is said to be 
asymptotically quantum ergodic if and only if 

lim Urn ^-^ '—- = V^ (64) 

for all acceptable sequences of observables {V/);eN- The limit Z — > oo plays the role of 
the semiclassical limit for quantum graphs. 

For the sequences of graphs satisfying ( |64] |. the rate of convergence is also of par- 
ticular interest. Therefore, we will treat a single finite quantum graph first, and come 
back to convergence and rate considerations afterwards. 

A calculation similar to (|6T1 ) shows that, for an observable V on G constant on each 
bond, one has 



(^\n"')= a 



sin(^vL) 
VL\\ + — -^— 

17 



a^)^{a'\VLW')+0{k-\ (65) 



1 V <^^^^ rfis^ 



where |a'') e y[ is the vector of ampHtudes defining ^P" through the construction ( fT2] l. 
There is a slight abuse of notation in this expression. On the left-hand side, V - 
(Vfc)fceNB is an observable constant on each bond, whereas on the right-hand side V 
stands for the diagonal 2B x 2B matrix Vbd.b'd' = SbM'^dM'^t- Such a matrix is called 
observable on yi and has mean value 

-_tr(yL) Zif 1 Vflh Zti VbLt ,_ 

This expression coincides with the mean value ( |59l ) of y seen as an observable on G 
constant on each bond. 

From dSOl l, (|65] | and (|66] |, we deduce that a quantum graph is quantum ergodic if 
and only if there exists a subsequence / 1-> v(/) of density 1 such that 

<VL),(,.) <a-<'>|VL|a''<')) 
hm s lim -:- -T^ = V (67) 

for any observable V on J{.. As above, one can restrict attention to observables V such 
that V -Q without loss of generality. 

A standard theorem of ergodic theory, proven for example in ||45ll , states that the 
quantum ergodicity property (|67] i is equivalent to the vanishing of 

Fv = lim ^ : 

for all observables V on yi with V = 0. Moreover, since the bond lengths are bounded 
by Lmin and Lmux by assumption (|62] |. this property is also equivalent to the vanishing 
of the fluctuations 

— ^ (VL)^(VLy,C,,, (69) 

for all observables V on yi with V = 0, where the intensity correlation matrix CjSjj' in 
the right-hand side is defined in i44i . 

In the case of an increasing sequence of graphs {(G;,5/));eN, asymptotic quantum 
ergodicity is obeyed if and only if the sequence {/^/^vJ/en whose terms are defined 
as in (l68l l, or equivalently the sequence {'Fi,Vi]ieN whose terms are defined as in i6% . 
converges to zero as / — » oo for all acceptable sequences of observables {V;)/eN- The 
rate of convergence is then called the rate of quantum ergodicity. 

The Gaussian Random Wave Models ( |48] l and ( |53] l predict the fluctuations 

as can easily be shown from the Gaussian predictions for the intensity correlation ma- 
trix ( |49] l and (|54] |. The parameter k was defined in ( l28T l. The term proportional to k 
describes the deviation from quantum ergodicity. For any admissible observable and 
bond lengths bounded by i6% the deviation predicted by the Gaussian Random Wave 
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Models is 0{B^^). Hence the Gaussian Random Wave Models predict that any increas- 
ing sequence of quantum graphs is asymptotically quantum ergodic and that the rate of 
convergence is larger by a factor of two if time-reversal symmetry is conserved. 
Note, that quantum ergodicity holds on average, in the sense that 



1 y <Vi)v 



Av^lim^^^2^^j^^V (71) 

ky<K 



for all observables V. This is known as the local Weyl law. It is easily checked to hold 
for any quantum graph . Indeed, by the definition ( l42t . Ay can also be written 



2B 



trL 



2B 



Yj<yi^)p{\aii\^)- (72) 



Then, the identity dHJ shows that /|a/jp) = {2BY\ and the definition (|66ll of V con- 
cludes the proof of the claim. The restriction to incommensurable bond lengths is not 
necessary for the local Weyl law, indeed we will show in the following subsection that 
W\\ is true for any choice of bond lengths. 

3.4 Green Matrices and Trace Formulae 

For (G, 5 ) a quantum graph, and for e > 0, one defines a sub-unitary quantum map 
Ueik) by 

U,(k) = T{k)S,T(k), with S, = e-'S, (73) 

and where T(k) is the propagation matrix of G given in ( fTSI l. The retarded Green matrix 
(resolvent) G{k) is the matrix-valued function on R+ defined by 

C«.(l-^,«r=^^^|S^. (74) 

It has poles in the lower complex half-plane at <p„{k) = Inp - ie for any p e Z. The 
advanced Green matrix G ' (A;) is the hermitian conjugate of G{k), that is 

Ct« = (,- „!«)-'= 2 T^J§^^ (75, 

It has poles in the upper complex half -plane at (p„{k) - 2np + ie for any p e Z. Making 
use of formula (IZSl l. it is not diflicult to check that, for any integer q > 2, and for any 
permutation cr e S q, the statistical quantities defined in dJTJ i with p = q read 

where, in the right-hand side, the average over k is defined by the formula 



1 r^ 

= Um - 

K^^ K Jo 



{f{k)), = lim - f(k)dk (77) 
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which is meaningful for any function / integrable on every compact interval [0,^]. 
The formula ( l76b relies on the non-degeneracy of the spectrum, which generically fol- 
lows from the incommensurability of the bond lengths. However, it still holds if the 
subsequence of levels ky that are degenerate is of density zero. There are other versions 
of the equality ( |76] ) where the right-hand side involve «, e N elements of G{k) and 
«a e N elements of G' (^) with «r + «« - q. A formula similar to (|76] | is used in Il20l 
to study the statistical properties of the eigenfunctions in disordered systems. For the 
derivation of exact expressions the choice of the permutation cr e S gin ( 1761 ) is mainly a 
matter of computational ease (and sometimes taste). Throughout the remainder of this 
subsection we will show how the different choices lead to different exact expressions. 
The Green matrices G(k) and G' (A:) can be viewed as the results of summing ge- 
ometrical series in Udk) and Ul(k). This gives rise to interpretations of their com- 
ponents as sums of walks on the quantum graph (G, 5). An oriented walk/? is a list 
(Po,/3i, . . . ,/3„) of consecutive directed bonds on the graph. Its topological length \p\ is 
the number of vertices traversed, that is \/3\ = «. The set of all oriented walks having 
topological length n is written W„. The metric length of yS is 

l(P)-^-^"f,L,^^. (78) 

The origin and terminus of y6 are respectively o/3 = Pq and tfi s p„. The set of walks in 
W„ having origin yS and terminus yS' is written W„(J3,/3'), and U„eNo W'„(j6,j6') = W(/3,/3'). 
We also define the stability amplitude 

n-i 

/=o 
With these definitions, it is easy to see that 

G(%,,= J] .-^\'*'<^U^ (80) 

peWlfi',13) 

and 

G\k),p.^ 2 ^'''^^'""%'- (81) 

Together with (|76] |, these formulae enable one to express the autocorrelation functions 
C[a] in ff2l i as sums over oriented walks. 

The diff'erent choices for the order of the left indices yS in (l76T l lead to different equiv- 
alent expressions for the autocorrelation functions Cia] in terms of oriented walks. In 
general, showing the equivalence between these trace formulae at the level of oriented 
walks turns out to be a very difficult problem. In this subsection, these non-trivial 
equivalences are illustrated by two alternative proofs of the local Weyl law (fTTT i. 

In the case of the intensity correlation matrix C^^', two permutations cr e ^2 of the 
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G 



^/3' 



E 



P' eW(f)',i3') 

1(0) = 0) 




E 



/?'£ VF(/3',/3) 
m = 10') 




Figure 2: The two equivalent formulae (l83T l and dSSl ) for the autocorrelation function 
Cjjij'. The underlying graph has not been represented for sake of clarity. The trace 
formulae (|88] ) and ( [89l ) are obtained from the ones represented here by adding the 
contributions where S and S ' are swapped and by dividing by two. 



left indices in (|76] | can be chosen. The identity permutation cr = id leads to 

Cpp' = lim|(G(^VG'(fcW)^ 

while choosing the transposition cr - (l 2) leads to 
C00' - lim|(G(^WG'(^V)^ 

= lim- y y e-(^>^'l>5,,.,,,.,/.A.,*. 



(82) 
(83) 

(84) 
(85) 



l3eW(l3,l3')P'eW(l3,l3'} 



In both cases, the Kronecker symbols originate from the average over k. 

These orbit expressions can also be recovered by means of the Poisson summation 
formula. If 6^{x) denotes the Lorentzian of width e centered at the origin, this formula 
leads to 



2B 



71=1 p=Q 

1 °° 



E? 



(86) 
(87) 



?=i 



The trace formula (|83T l, or more exactly its symmetrization obtained by replacing 
AaAa* with jiAgA^* + A/Aa), follows from ( |87] i and the identity 



Cflfl' 



27rV 
B 



{@e{k)l3l3@e{k)l3'l3')^ 



(88) 
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This identity is a consequence of the fact that, in terms of distributions, the product 
2n£6^{x)6eiy) tends to zero if x ^ y and to 6(x) if jc - y. Similarly, the symmetric 
version of dSST l follows from dSTl ) and 

Cpp' = ^ {0e(k)p'p 0e{k)pp)^ . (89) 



The main advantage of the expressions ( [88] l and (|89] l, involving Qe{k), over their ana- 
logues ( |82] | and ( l84l i. which involve G(A:), is that the matrix ^^C^) is real, whereas G{k) 
and G' (A:) have non-vanishing imaginary parts and must always appear together in ( |76] |. 
In particular, the first moment M^j - (|a/3p ) can be written 

Mp,,^^\im{g,{k)pp)^^ (90) 

which involves a single closed oriented walk, while, in terms of matrices G{k), an 
additional directed bondjS' must first be introduced in order for Mp^i to be written as 

the sum Hb^i Cpp' - Ym=\ (|fl;Sp|a;S'P) ^nd the representations (|83] ) or dSSb to be used. 
From ( |87] ) and ( |90l ). one finds directly that 

%,i = ^, (91) 

which, together with ( l72l i. provides a second proof of the local Weyl law. Let us now 
use the trace formula dSSl l and perform the sum over the directed bondyS'. It is easy to 
show by induction over « and m that the unitarity of the scattering matrix S implies 

IB 

Z Zi Yj ^ 10110' Ai^^0'* ^^"''" (92) 



for all n, m e Nq. With (l85l l, this gives 

2B oo 

M,, = J] C,p, = lim I 2 .-<-'">5„,„, = -, (93) 

/J' = l «,m=0 

which, together with ( |72] |. yields a third proof of the local Weyl law. 

The choice for the permutation cr e Sgin ( |76] | leads, in the case q = 2,to the equiv- 
alent expressions ( |83] | and dSSl l for C/j/j- in terms of oriented walks that are illustrated in 
Figure|2] Similar pictures could also be drawn for C|^„ . .^^ ,] when q > 2 . Indeed, the 
right-hand side of (|76] i with B'. - B: for all < / < a - 1 and with a fixed permutation 

cr e S q can be expressed as a sum over ^ oriented walks ySi , . . .,Pq, where each ySy leads 
from the directed bondySy to the directed bondy6o-(y). The waUcp^ is followed with Sl, 
whereas the q - I other walks are followed with S f, and its metric length must equal 
the sum of the metric lengths of the q - I other walks. Therefore, (|76] | yields q ■ q\ 
different ways of expressing the autocorrelation function C[^] of degree q in terms of 
oriented walks. Here, the first factor q accounts for the q possible choices for the walk 
followed with S I . 
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3.5 Long Diagonal Orbits 

An expression for the fluctuations Ty in ( |69] l can be obtained by retaining only a subset 
of the whole set of pairs of oriented walks 0,p') entering ( [83T l. For this purpose, it is 
convenient to come back to the expression dSSl l of Cpp' in terms of ^^(A:) and write the 
fluctuations of an observable V with V = as 

The right-hand side can be written in terms of periodic orbits rather than closed oriented 
walks as in ( l83b . A periodic orbit is an equivalence class of closed oriented walks 
whose sequences of directed bonds diff'er from each other by cyclic permutations. For 
a periodic orbit p, the notions of reverse /3, topological length \p\, metric length Ip 
and stability amplitude Ap are inherited from the oriented walks terminology, and the 
repetition number r^ is the number of times p retraces itself. With this notation, one 
gets from dSTJ l 

tr {0,VL) ^-%y e-l/'l iXf:^e'*'.A„, (95) 

where the sum is over all the periodic orbits on the graph and {VL)p stands for the 
number obtained by accumulating the values (VL)/? of VL along p. The square of the 
last formula admits the spectral average 

([tr(^.yL)]^) =-L y ^ZfMXf:)i5l(A,A;).-<l"l-M). (96) 

p.q:li,=l,i ' ^ 

The diagonal approximation, which consists in only keeping the pairs q = p and q = p 
in the time-reversal invariant case, yields 

([tr(^.yL)]^r = ^ X ^1^.1 V^^'^ (97) 

p,q:lp=l,i P 

where k is the parameter as in ( |70] i indicating whether time-reversal invariance is bro- 
ken or conserved. We have neglected some corrections in the diagonal approximation 
which are due to repetitions and self-retracing orbits. These can be shown not to con- 
tribute in the present context. The formula ( |97] l is then approximated further The 
orbits for which r,, > 1 are rare, so that we only keep the primitive orbits, namely 
those with rp - \. We also take the long orbits approximation Iil8i , which amounts to 
approximating 

[{VDpf ^ [{VL)X ^ \P\'^^. (98) 



Besides, the stability ampUtude is known to behave like ||25]| 

|A„|2 ~ e-«l'^l, (99) 
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where a is the topological entropy. This parameter also characterizes the number 
|^|-iga|pl of periodic orbits having topological length \p\. With all these approxima- 
tions, (|97T i reduces to the integral 



([^(^^-)nr^iX" 



^"''d\p\ ^,^^tr(VLf ^^_,|, 



-2e\p\ 



\P\ 



2B 



Hence 



r: 



diag^ trjVLf 



(trLY 



(100) 



(101) 



This formula, obtained from the long diagonal orbits, coincides with the prediction 
dTOl i of the Gaussian Random Wave Models ( |48] | and ( |53] |. It predicts asymptotic quan- 
tum ergodicity for any increasing sequence of quantum graphs and a universal rate of 
convergence B\ as in IfTSl . 



4 Generating Functions 
4.1 Definition and Principles 

The Green matrices introduced in the subsection |3.4| can be obtained as the derivatives 
of certain determinants. It is convenient first to introduce a Grassmann algebra A, 
which can be decomposed as the direct sum of its commuting sub-algebra Ag, called 
bosonic, and a set Af of elements anticommuting with each other, called fermionic. 
Then, the amplitude space ^ can be graded to get ^©^, and the Grassmann envelope 
(J?l ffi .^)(A) defined as in fSl can be built. This set reads 



(^ © ^)(A) = \y = iyl]-^yB/F^Yj ^i/ph^^ K/F e ^B/F 



(102) 



where the elements |eyj) refer to the elements in ( fTTT i of the natural basis of yi. The 
elements of (^ © ^)(A) are called supervectors. The set of endomorphisms on (I102l i. 
once written in the natural basis of yi, form a set of supermatrices written L{J\\J[). For 
g' > 2 an integer, let us introduce complex numbers y'l, . . . , yV/-i and jq, respectively 
referred to as retarded and advanced sources, and let us also consider q directed bonds 
a\, . . . Uq-x and ao- The corresponding retarded and advanced source supermatrices are 
defined by 



JrUr) = l+EB^jrE^'K 
JaUa) = l+EB®jaE'-"\ 

where Eb is the projector onto the bosonic sector of (^ © ^){A), 



(103) 
(104) 



J = 



Ja 
jr 



JO 
j\ 

^ jg-l ) 






£00. ao \^ 
£ai.ai 



(105) 
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and, for any two directed bonds a, a' € N2B, E"'"' stands for the 2B x 2B matrix whose 
components are {£"'■" )pi}' = 6a.pSa'.i3' in the natural basis of J{.. The number q - I 
of retarded sources corresponds to the number of matrices E"''"' contained in £'*''', 
so that the product in ( |lQ3t makes sense. In (11031 ). ( |104t and in what follows, some 
unit matrices or supermatrices are not explicitly written in order to keep the notation 
as simple as possible. For example, the symbols 1 in ( 11031 ) and ( I104l i must be read 
\BF®^y{, where Ibf is the unit supermatrix in Bose-Fermi space and Ij^ is the 2Bx2B 
unit matrix in amplitude space ^. 

Let q > 2 and let [or] = {a{),ai, . . . ,aq-\\ be a list of q directed bonds. The 
corresponding generating function is defined by 

f[„]C/) = (sdet-l(l - J,Ur) ■ U,{k))[\ - Jaija) ' Ul{k))) ^ , (106) 

where y^C/V) and Jaija) are defined fr omj s {jajrf = (;'o, ;'i, ■ ■ ■ , jq-iY and from the 
directed bonds in [a\ as in ( |103t and ( |104| i. Notice that this function is well defined in 
a neighborhood of the origin, and that it also reads 

^[„]C/) = (det-'(l -jrE^'\G{k) - l))(l - uE^"\G\k) - 1)))^ (107) 

in terms of Green matrices. 

It is convenient at this point to give a general rule governing derivatives of determi- 
nants of the form ( 1107b . An important quantity is the p factor 

/ \ number of cycles in o" /iao\ 

Pa((r) = a ^ (108) 

defined for any a e K and any permutation cr e S'., of i e N elements. This factor 
can be seen as a generalization of the signature (-1)°" of cr e S s since the identity 
(-1)°" = (-l)*p_i(cr) holds. Now, if A = (A*", . . . , A'")^ is a vector containing 5 e N 
square matrices A*'' of size « g N and if 7 e C*, we have the equality 



dji . . . dj. 



detd-jA)- 



>-0 ireS, /=1 Xi=l 



This result can be proved by induction over s. The right-hand side has a natural dia- 
grammatic representation where each / e N, is a point and where an arrow is drawn 
from i to j whenever cr(/) - j. The sum in ( |109t is then the sum over all such diagrams 
in which each point / € N.s has exactly one outgoing and one incoming arrow. The 
value of each diagram is a product of traces of the type tr (A*'*A°"''' ■ ■ ■ A°"''*'^j, with p 
being the smallest number in No such that cr''+'(/) - /, weighted by its p factor, which 
can be deduced from the number of connected sub-diagrams. 

Let q > 2 and let [a] = [ao, ai, . . ., a^-i] be a list of q directed bonds. The rule 
( |109t can be applied to the expression ( |107| l for the generating function, and, making 
use of (l76T l, one easily gets 
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where Cia] is the autocorrelation function defined in ff2] l. and 



^^[ff] = 






^[a] 



(111) 

7=0 



The denominator (^- 1)! in jllOb comes from the number of diagrams arising when the 
rule (1109b is applied to the q-l retarded derivatives on ^[a](j)- By (ESI), these diagrams 
all yield the same contribution. 

It is not difficult to check that the generating functions have the following property. 
For all ja and jr - {j\, ■ ■ ■, jq-i) in a sufficiently small neighborhood of the origin, 

f[«]Oa,0)=^M(0,jV)=l. (112) 

For cr e S g and [a] = [ao, ai, . . ., a^-i], one can introduce a function ^[^.(j) by the 
formula ( 1106b using the matrices E"''"'^'^ in place of £■"''"' in the source supermatrices, 
and (IllOb then serves as a definition for Cf,. The function ^[^,(j) also satisfies the 
property ( 11121 ). and the identities ( |76] | and ( |109b ensure that Cj^, = C[a] for any cr e S q. 
In what follows, the arbitrary choice for cr e S q va ^f^,(j) will be called the choice 
of convention. These different but equivalent expressions must not be confused with 
the equivalent sums over orientated walks which are the object of Subsection l3.4l Any 
convention cr e 5 ^ for the generating function involves {q - 1 ) ! equivalent sums over 
orientated walks. However, in the case q - 2, the permutations cr = id and cr - {0 1), 
which are referred to as parallel and crossed conventions in the sequel, do correspond 
to the sums ( 1831 ) and dSSl ) respectively. 

We started this chapter with the convention to choose o" e 5^^ to be the identity and 
we will use this convention in most of the following calculations. This convention is not 
only singled out by simplicity; it results in a generating function ( 11061 ) that is explicitly 
gauge invariant while in other choices the gauge invariance is only restored in the limit 
dllOb . It also reduces the complexity of some calculations because the matrices JAJr) 
and Jaija) for the source terms ( 11031 ) are diagonal matrices. While each convention 
yields a different but exactly equivalent expression approximation schemes may break 
the exact identity. This is not worrying as long as the difference is in sub-leading order. 
For time-reversal invariant graphs the generating function ( 1106b is usually not explicitly 
invariant when one replaces any directed bond by its reversed partner for q > 3. The 
invariance is only revealed once the derivative in ( II 10b is taken (the limit e ^ is not 
required). 

4.2 Diagonal Approximation 

Before working further on the generating function ( 1106b with the supersymmetry method 
we introduce in this subsection similar generating functions and develop a correspond- 
ing trace formula. The diagonal approximation to this new type of generating functions 
turns out to behave very differently from the oriented walk representations previously 
discussed in the subsections l3.4l and l3.5l 

The definition ( 11061 ) of the generating functions, and the fundamental formula dllOb 
can easily be generalized to any correlation function dJTJ i with p - q. Moreover, these 
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correlation functions can also be written in terms of logarithmic derivatives with some 
analogy to ( II 101 ). We will focus on the case p - q -2 for which the general correlation 
function can be written as 

{a*a^aa'^*aa,aa'^j = lim — 6E[a^,a\;a2,a'^], (113) 

where 

^la„a[;a„a',]Ua,jr) = (logdet(l - AOV)^/.^) log dct (l - JaUaWl(k)))^ 

= ( log det (l - lUrWAk)) log det (l - /,(;J^C/,(A:))')^,(1 14) 

the source terms are given by 

lUr) ^ I + jrE"''"'- and Uj,) ^ 1 + j„E">K (115) 

and 6E = ^"g. E.(ja,jr)\ ■ The intensity correlation matrix can be obtained in two 

different ways 

Caa' = ^i^€^O^SE[a,a':a,a'] (116) 

= \im,^0^6Eia,a':a',a] (117) 

referred as parallel and crossed conventions, respectively. In the orthogonal class one 
has a third representation 

Cac' = lim — (5H[q,_q,;q,',q,'] (118) 

called time-reversed crossed convention. 

The formula log det = tr log enables us to write the new generating function jlMI l 
in terms of generaUzed periodic orbits on the graph. Indeed, expanding the logarithms 
and performing the spectral average yields the trace formula 



°° 1 

^lai,a'^;a2,a'^]UaJr)^ ^ ^ Yj ^"^npUrf [Aa.p'UaTf Spl^^^p'l^, (119) 

where the retarded and advanced modified stabihty ampUtudes Ar,p{jr) and Aa,p(ja) of 
the generalized periodic orbit p - /3iP2- ■ ■P\p\ are defined by 

\P\ \P\ 

A,..p{jr)=Y[{lU,)Se\p,^,p, and Aa.p{Ja) = Y[V''(JafS,]p_^^p_. (120) 

The periodic orbits p and p' in ( I119l l are all primitive but can be of a slightly more 
general type than the primitive periodic orbits considered in Subsection 13.51 Indeed, 
the retarded source term in ( 1115b introduces the possibility to jump 02 r\ ai, and 
similarly the advanced source term introduces the possibility to jump o-'j rv ai,. The 
set Po,a2 in ( 11191 ) then contains all the primitive periodic orbits that are compatible 
with the topology of the graph with an additional bridge a2 r\ ai at the center of these 
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directed bonds. Note that the two sets Pq.,,,, and P„;„' of generalized periodic orbits 
need not be identical. In the parallel convention for Caa' , the source terms are diagonal 
matrices, and Paa = P reduces to the set of standard primitive periodic orbits, which 
only respect the topology of the graph. The length of a generalized periodic orbit is 
just the sum of all bond lengths along the periodic orbit, where every jump a2 r\ o-i 
contributes ^(Lq., + La,). In the trace formula ( |119l l, only pairs of primitive orbits 
contribute such that a repetition of one orbit has the same length as a repetition of the 
other 

The diagonal approximation to the trace formula ( I119l l reduces the sum over pairs 
of primitive orbits to either equal orbits p' = p, or time-reversed orbits p' = p. In both 
cases, the factor Spi ^p'l^, enforces p - p'. Note that p has to be a periodic orbit in the 
intersection p e Paia2 ^ Pa'^o' to contribute to the diagonal approximation. 

The remaining sum over periodic orbits can be resummed 

^^diag , . . -. ^diag.D , . . , .-^diag,C / • • \ /^ 

- - log det (l - M^'ija, jr)) - log det (l - M'^Ua, jr)) + C 

(121) 

where M'^{ja,jr) and M'-{ja,jr) are modifications of the classical map ( l32T i. They 
describe diffiison and cooperon propagations, which originate from pairs of periodic 
orbits with p' - p and p' - p, and are given by 

M (ja,ji-)l3,l32 = / J J r(ji)pili' Ja(ja)l3'l3i\S e. Ii'li2\ 
P' 

= (1 + jrSl3ia,Saia2 + ja^/Jia'^Sa'y^j \S e. /3,/32 I + 

jrja^fi, a, Sa,ai, Sa\ o-j 1'^ f . aj/S: I ( 1 22) 

M {ja,jr)l3il32 = / , Jr{jr)fiil3' Ja{ja)i3'i},S ^_ p'p, [S ^ g,^J 
P' 

jrja5|3^cn5a^&[6a2&\S e, a2P2\^ c,ii2a2) ' (123) 

The term C in ( 11211 ) contains corrections for repetitions and self-retracing orbits {p - p) 
which can be shown not to contribute to our final result and will be omitted henceforth. 
Recall that the classical map is defined by Mp^|3, - |5^,^, p, so that M^(0, 0) - M^ = 
e^^'^M. For time-reversal invariant systems, the products of scattering matrices in ( I123l l 
reducesto5f,/s,/j2(5^j,^_)* = \Se,p,p2\^ = ■^f./SiA- so that limf^o ^^(0, 0) = M. If time 
reversal symmetry is broken M''(0, 0) does not reduce to M and it does not describe a 
Markof process on the graph. We will see that the cooperon term only contributes to 
the diagonal approximation formula if time-reversal symmetry holds. 



28 



The derivatives with respect to yV and ja can now be taken and yield 



CI— idiag,D 



tr 



1 d^M'^ \ I 1 dM" 1 dM'^ 
+ tr 



~^a]a'-,^a',t 



1 - M, dudj 



1 - M, dj^ 1 - M, dj- 



j.,=jr=0 



l-M, 






(124) 



for the difFuson generating function, and 






M^(0,0) 



M'^(0,0) \ / M^{0,0) 



1-MC(0,0)L ,,\1-MC(0,0)L 



(125) 



for the cooperon generating function. For broken time-reversal invariance the classi- 
cal cooperon map M*- (0, 0) has no unit eigenvalues in the hmit e — > 0, and hence, 
(I125l l identically vanishes in that limit. By contrast, in time-reversal invariant systems 
M'' (0, 0) = Mf and the cooperon generating function does contribute. 

Finally, only terms in ( 11241 ) and ( 11251 ) that are singular as e — > contribute to 
the correlation function (I113l l. In order to isolate these terms, one makes use of the 
decomposition (l35T l of classical orbits as the sum of a uniform component |1)(1| and a 
massive part R. This yields 



_ ;^(l-2e) 



-5e: 



Oa]a2^a''^a'-, ' . ^9 \^<y\aij^Q!'a2 ' \^ ^)^aia'^a2ai,j' 



+ J^Sa^a.Sa^a'^ {[Ra,a\ + Ra\a,] + (^ - 1) [Ra^a] + R&^a^]) + 0(e) 

(126) 



After dropping terms that are 0(e), ( 11261 ) may be expected to provide an approximation 
to the generating function ( 11131 ). However, its first term diverges like e"' as e ^ 0. 
At first sight, this seems to make periodic -orbit analysis using the trace formula ( |119t 
for the generating function much less useful than the previous trace formulae from 
Section 13.41 which behave nicely in the diagonal approximation. On the other hand, 
the same divergence also occurs in the analysis of spectral correlations, which become 
singular in the diagonal approximation at small energy differences. Indeed, one may 
obtain the corresponding trace formula for the spectral two point correlation function 
on a graph R2(s) by replacing the source terms Ja(ja) and JAjr) appropriately. In this 
context, a supersymmetry method developed in Il22l|23l . which will be adapted to our 
purposes in what follows, cures the divergence. One may also try to add off-diagonal 
terms in the trace formula in a systematic way but we will not pursue this here. 

Note that the Kronecker symbols in the expression ( 1126b force the correlation func- 
tion (a* a* iQa.aa' ) to vanish for all combinations that are not invariant under all local 
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gauge transformations allowed by the unitary or orthogonal symmetry class. The non- 
vanishing combinations are then equivalent to the three different conventions (1116b . 
([TT7J and (III8I 1 for expressing the intensity correlation matrix. These three conven- 
tions lead however to different formulae. 

For any of the three conventions ( 11161 ), ( |117| i and dllSI ), if one replaces e •-» ^, the 
first line in ( 1126b reproduces the prediction of the Gaussian Random Wave Model up to 
corrections that are 0{B^^). The second line in ( 11261 ) then gives a correction in terms 
of system dependent massive modes. This massive correction turns out to be different 
for the three conventions of expressing Caa' ■ 

In the orthogonal class only the parallel convention ( 1116b provides an approximated 
intensity correlation matrix that respects the identity Caa' - Caa' satisfied by the ex- 
act intensity correlation matrix. By contrast, if either the crossed convention ( 1117b 
or the time-reversed crossed convention ( 11181 ) is used, the massive terms in the ap- 
proximated intensity correlation matrix explicitly violate this symmetry. The origin of 
this discrepancy is that, with the parallel convention (11161 ), each of the two logarithms 
in the generating function 'B.[a^a';a,a']{ia,ir) is invariant under time inversion, while in 
the crossed and time-reversed crossed conventions the symmetry is only restored after 
taking the derivatives and performing the limit e ^ 0. 

The observations above concerning time-reversal symmetry makes the parallel con- 
vention ( 1116b a privileged choice when it comes to the diagonal approximation. This 
convention yields 

Cf J'" = i 1 ''^^"^^^ + Saa' +(k- 1)6 aa' + Raa' + Ra'a + (« ' I) (Raa' + R&'a) 



4B2 \ 4B£ 

(127) 
The three first terms are universal, and are equal to the prediction of the Gaussian Ran- 
dom Wave Model if e is chosen finite and set equal to ^ (which we cannot justify at 
this stage). The remaining three terms involve the matrix R and they describe massive 
corrections to the universal result. In fact, these massive contributions may dominate 
the correlation functions and, as a consequence, the rate of convergence for quantum 
ergodicity, or they may destroy quantum ergodicity altogether. This point will be dis- 
cussed further in Section|2l 

4.3 Nonlinear Supersymmetric cr Model 

The generating functions in ( |106b depend strongly on whether time-reversal symmetry 
is broken or conserved. Time inversion acts on supervectors if/ in the Grassmann en- 
velope {X e A')(A), defined from A" = :?[ ® C" for some n e No as in ( |102b , and on 
supermatrices A € L{X\X) as 

Tif/ = crfi/r* and A'^ = o-fA^cr'/. (128) 

In ( 1128b , CTj' is the first Pauli matrix acting on the direction space ^j, tfr* denotes the 
vector obtained from tfr by taking the complex conjugates of each component, and A^ 
is the transpose of A defined as in |fT9l by the condition {AijJiY\jj2 - il/\A^ijj2 for all 
^Ai , ^^2 in (/Y © A')(A). Here and henceforth, if/ stands for the row vector obtained from 
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the column vector ip e {X ® X){A) by usual transposition. One can now introduce a 
2-dimensional C-linear space TR, the time-reversal space, and the mapping 

from {X®X)(A) to (XiS>TR®XiS>TR){A) called time-reversal doubUng. We work with 
the convention;^** = -x for all;i' e A/r as in |fT9l . and hence, the Hermitian conjugate 
of ^P in dug reads 

T = («A'",«AVfcrf), (130) 

where i/'' = i//*^ is the Hermitian conjugate of ifr, and cr*^ stands for the third Pauli 
matrix acting on the Bose-Fermi space. Similarly, the time-reversal doubling of a su- 
permatrix A e L{X\X) is defined by 

^"^^(« ;o„=(o <.v,')„ 

and is an element of L{X ® TR\X ® TR). In ( |129t . ( 1131b and in what follows, an index 
TR added to a supermatrix means that this supermatrix is explicitly written in the TR 
space, and the same notational trick is used for any other space. The components in 
time-reversal space will be indexed by f e (t, i). The definitions above call for a notion 
of generalized transposition yC of yi e L{X ® TR\X ® TR), which is defined as in ||23l 
by 

/ A BF \ 

y[^ST^V, where t s erf , ^ ■ (132) 

This definition implies that the equality *i'iy[*I'2 = *i'2^^*I'i holds for any couple of 
supervectors ^Pi, ^2 e (A" » T/? e A" » TR)(A) and for any supermatrix Jl 6 L{X ® 
TR\X(S>TR). It follows that (Jl^ = ^ and (^IS)'" = S^Jl" for any such supermatrices. 
Moreover, using the property (A^)^ - cr^^Acr^^ of the transposition in L{X\X), it is 
easy to check that a supermatrix J?l obtained from some A e L{X\X) by time-reversal 
doubling ( I131l l is invariant under generalized transposition. 
Now, the generating functions can be written 

^mC/) = sdet-'7,y„ (sdet-i(j',:i - 'U,{k)){j;' - Vlik)))^ , (133) 

where J'rja and lAeik) are the time-reversal doubles of 7,./a and U^{k). Following the 
scheme developed in Il22l and Il23l . the generating functions ( 1133b can be represented 
in terms of a nonlinear supersymmetric cr model. 
First, it is convenient to make use of the equality 

A B \ , . , .., .-lDr.-l/ 



sdet M; " =sdet(A/))sdet(l-A-'BD-'C) (134) 

that holds for any square supermatrices A,B,C and D of the same size, and write the 
retarded and advanced superdeterminants in ( |133b as 

sdet-'/2(j-;i _ .j/^) ^ sdet-'/2 ( y, ^ ^,^ j , (135) 
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and 



sder'/^(j-;' - 'Z/J) = sder>/2 1 . ^^„ ^'7' 1 • (136) 



r^' V3:'- J- 



In these expressions, S^ is the time-reversal double of S e, and the square root of a 
matrix is defined by keeping the same eigenvectors and by taking the square roots of the 
eigenvalues fixing the half-line singularity of the logarithm to (-00, 0]. It is not difficult 
to check that this definition of the square root leads to the natural properties a/A VA - A 
and yA^ - Va' . These two properties have been used in order to obtain ( 1135b and 
( |136l l. Besides, if ^ is the time-reversal double of A, then V^ is the time-reversal 
double of Va. The 2-dimensional structure introduced in the right-hand sides of ( I135l l 
and (1136b is referred to as the auxiliary space X, and the components with respect to 
the basis of X used in these two matrix expressions will be indexed by x e (1,2). The 
formulae (1133b . ( 11351 ) and ( 1136b now enable us to express the generating functions as 
the Gaussian superintegrals 

f[„]C/) = sdet-'7,y, ■ fd^{e-'^''^)^ , (137) 

where d^) - dijjrdif/a, dif/r and d\jja being two Berezin measures fS] on (^ © y[)(A), 'V 
is the time-reversal double of ijj, and 



.™M*..^.kfn(^;;) 



■(*.*.)(,, u "fTil-iy (.38, 



The indices r and a of iff and *P refer to the retarded and advanced components of these 
supervectors used to write ( 11351 ) and ( 1136b as Gaussian superintegrals respectively. 

In the quadratic form ( 11381 ), the off-diagonal couplings depend on k through the 
variables (kLi, . . . , kLg) in the propagation matrix T. Using the fact that T and yfS^ 
are invariant under generalized transposition, these off-diagonal terms can be written 

5ef[T] = 2>I',i ^/s',T^,2 + 2^a2T'^- ^/s',^'^'al . (139) 

Since the bond lengths are assumed incommensurate the invariant measure of the auto- 
morphism k i-> (kLi , . . . , kLg) (mod 2n) on the B-torus is merely the product of B Haar 
measures on the circle l?). Hence, 

(e-^-[*])^ = f] r ^e-s^-rm^**], (140) 

where 

S'fVi'b-, 'Pb] = 2Y, [(%■! V^),„, e''^"^',.2M + ^aiMde-'"" ( V^"'"^«4J . (141) 
d=± 

Then, the color-flavor transformation lISTl can be applied separately to each integral in 
the right-hand side of (|140b . This procedure introduces supermatrix variables Z/, and 
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Zb, which lie in IX^d ® TR\^d ® TVi), and yields 



with 



(,-^c.[>r,^^ ^ p] r^(z,,z,)sdet(l - Z,Z,).-^crt*'^Z"-^"I, (142) 

5^f[>Pi; Zi, Zi] = 2 2 (Th VisT)^^^ Zfc;rfrf- ( 4Se^^a\^^, 

d=± 
+2Y,^a2-bdZb,dd''i'r2-bd (143) 



Basically, the retarded and advanced components of 'I', which are uncoupled in ( I141I I 
become coupled in ( 11431 ). and conversely, the components in auxiliary space, which 
are mixed in (11411 ). are diagonalized by the color-flavor transformation. The reason for 
resorting to this transformation is to get an action with saddle-points, which is not the 
case in (1141b . The integration in (1142b must be performed over the set of supermatrices 
{Zb,Zb) satisfying the conditions 

Zbb = Zl^, Zff = -zip, (144) 

and such that the eigenvalues of the positive hermitian matrix Z^^Zbe are less than 
unity. The measure d(Zb,Zb) is then the Berezin measure over this set. 

In order to simplify the notation, one can introduce the new supermatrix fields 

B B 

Z^^Zh and Z^QZb, (145) 

b=\ b=\ 

which belong to L{J[ ® TR\J[ ® TR). These supermatrices still satisfy the color-flavor 
requirements (1144b . From (1143b and the diagonal terms of (1138b . one gets the new 
quadratic form 

The integral over ^ in ( 1137b remains Gaussian after the color-flavor transformation, 
and, from the exphcit formula ( 1146b . the generating functions become 

^mC/) = sdet^'y,^ rt/(Z,Z)sdet(l -ZZ)sdet"'^2| ^ ^_A 

The first superdeterminant in the integrand comes from the B superdeterminant factors 
introduced in (1142b . Making use of the rule ( 1134b once again and resorting to the 
well-known formula sdet - exp str log enables us to write 

^la^U) ^ ^ d(Z,Z) e-'^^-^\ (148) 
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where the function S [Z, Z], called the action, or the exact action in order to distinguish 
between S [Z, Z] and its subsequent approximations, is defined by 

S[Z,Z] = -strlog(l -ZZ) + -siilogil-ZSlZ^S^) 

+ istrlog(l-J-,Z^J-„Z). (149) 

Notice that, if the sources jV and y„ are set to zero, the resulting source-free action 
So[Z,Z] is precisely the one obtained in Il22ll and 1231 for the generating function of 
the spectral two-point correlation function. 

The different conventions cr e S q for the generating functions discussed at the end 
of Subsection l4. 1 l ean also be written in terms of the nonlinear supersymmetric cr model 
(I148l l. Indeed, in order to get ^f^iC/)' it suffices to replace 7, and Ja with 7,°" and J'^ in 
the exact action, where these two new source supermatrices are defined as in ( 11031 ) and 
( |104t using the matrices E"''"''''' instead of E"''-"' for allQ < j < q - I. 



5 Mean Field Theory 
5.1 The Zero Mode 

The first step of our approximation scheme consists of restricting the superintegral 
( |148l l to the subset of supermatrices (Zq, Zq) around which the first variations 

5o[Zo + /7W,Z„]-5o[Z„,Zo] . ,. So[Zo,Zo + t]W]-So[Zo,Zq] 
hm and hm (150) 

;;^0 J] n^O J] 

of the exact source-free action So[Z,Z] vanish as e — > for all supermatrices W in 
L(TR ® ^\TR ® J?[). This subset of mean field configurations, called the zero mode, 
was identified in ll22l and ll23l and consists of the supermatrices satisfying 

Zo = Ij^ ® y and Zq = t:n ® Y, 
with Y,YeL{TR\TR) such that Y ^ Y\ (151) 

Moreover, Y and Y must be diagonal in TR space if time-reversal symmetry is broken. 
Of course, the color-flavor relations (1144b must still be satisfied, that is, the identities 
?BB - igB and ?ff - "^'pp are fulfilled, and the eigenvalues of IgglfiB must have mod- 
uU smaller than one. The supermatrices (Y, Y) satisfying these relations parametrize a 
manifold, the so-called Efetov cr model space. Efetov's cr model space with unitary 
symmetry has 4 commuting and 4 anticommuting parameters, whereas 8 commuting 
and 8 anticommuting parameters are involved in the orthogonal symmetry class. 

Let us introduce a 2-dimensional C-linear space RA, called retarded-advanced space, 
and let us consider the supermatrices in L{RA (8> TR\RA (8> TR) 

R^ll^] and /?-' = ( _T^ "^ ] , (152) 

\ Ira \ i-YY i-YY Ira 
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Then, one can set lISTI 

Q = Ro^'^R-\ (153) 

where cr^"* stands for the third Pauh matrix in retarded-advanced space. By construc- 
tion, these matrices satisfy Q~ = Q. Moreover, if for a supermatrix A having a retarded- 
advanced structure A denotes the supermatrix 

A='KA^'K, where '^ = ( q \\ ^ (154) 

the Efetov cr model space is characterized by the constraints Q - Q, Q^ - cr^Qcr^, 
and Q diagonal in TR space for the unitary symmetry class. Efetov's polar coordinates 
|fT9l then involve writing 

Q = UQqU, (155) 

with 

cos / sin \ A _ / i^B 



I -'Sine -COS0 /^^ \ Of j^^ 

The equations Q^ = 2o and go = Qq are automatically fulfilled for any real symmetric 
matrices 6b and 6f acting on the TR space if 

™(«"b(«:L(o:L -' 

are required to satisfy UiUi - 1 and U2U2 - 1, that is S,- = u- - uj^ and v, = 
cr^^v-cr^^ - v^', for / € {1,2). The purpose of Ui is to diagonalize Q in Bose-Fermi 
space, and hence, this supermatrix contains all the anticommuting parameters. One can 
for example choose 

with 77,, 77* € A/r, / € {1,2), and define vi by substituting iki for 771 and k2 for 772. For 
the Efetov space with unitary symmetry, one sets 772, a:2 — > 0, in which case 7777*77 and 
77' 7777' vanish, and similarly for k and tc' . Requiring go to carry the additional symmetry 
Qq ~ cr^'^Goo^'* amounts to writing the matrix angles 6b and 6f in (1156b as 

with 01 , 02 > and 6 € [0, 27r], and set 02 ^* in the unitary symmetry case. Together 
with the property C/[ = Ui which follows from the definitions of ui and vi above, 
and from imposing C/j = U2 on U2, this symmetry implies that the required equality 
Q^ - cr^'^gcr^'* indeed holds. There are still 2 and 5 remaining commuting parameters 
that have to be included in U2 in order to span the full Efetov space for unitary and 
orthogonal symmetries respectively. It is not difficult to check that, for any matrix V in 
SU(2), and for any ^,x e [0, In], 

«2 = n T, and V2 s \, : (160) 



^ I BF \ ^ I BF 
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lead to a supermatrix Ui in (I157l l with the required symmetries £72 = f^J ~ ^2 '■ ^^^^ 
finishes the description of the Efetov space with orthogonal symmetry in terms of polar 
coordinates. In the unitary symmetry case, one can take;^' = and V - e"^""i . 



5.2 The Mean Field Autocorrelation Functions 

The restriction of the superintegral ( 11481 ) onto the zero mode defines mean field gen- 
erating functions ^^Aj) which, using ( 11521 ) and ( 1153b . and after some algebra, can be 
put on the form 

&)^\dQe-'TiQ^P^^^^), (161) 

where 



Be 



5r[e] = -rStrG + ae) 



(162) 



is the source-free action ^o at the configuration Q of the zero mode, Q denotes the 
supermatrix cr^Q — 1, Q^ stands for its Bose-Bose block, and P[a]{j) is the supersym- 
metry breaking factor 



P^a]U) = det 



•4, 



£(r) 







Q' 



RA 



(163) 



In this last expression £'"* is the time-reversal double of £'"* in ( I105l l. and similarly 
£*''' is the vector containing the time-reversal doubles of the q-\ matrices entering the 
vector £''''', and in (1161b . dQ is the measure d{Z, Z) in ( 11481 ) induced on the zero mode 
manifold. Notice that the scattering matrix S does not enter the mean field generating 
function. It can indeed be seen in ( 1149b that, after being commuted with Zq and Zq, S 
meets its adjoint in the mean field action and thus disappears by unitarity. 

The formula ( II 101 ) applied to the mean field generating functions instead of the 
exact ones generates mean field autocorrelation functions Cf^p and commuting the 
derivatives with respect to the sources with the superintegral in ( 1161b yields 



C, 



'MF 



lim 



(2e)g-' r 
•B(q- 1)1 J 



^0 2B(q 



dQe'^o 6P[, 



(164) 



where 6P\n} denotes the derivatives 



SP[a-\ = 



(q-l „ \ 



*..!=0 



djs 



P\a0)- 



(165) 



These derivatives can easily be calculated by means of the general rule ( 1109b . For any 
integer q >2, one gets 



a-eS,, ' <elT,il' 



(166) 
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where, for a € K and cr e S q, Paicr) denotes the p factor defined in ( IIO8I 1, and for any 
vector t 6 {t, W, 

i Gf„;,„,,„n^;;e?,,,,.„ ifcr(o) = o 

"^''^^" Qiv..,..&,.,U^Q^n,,.,, ifcr(0 = 0,/6N,_, (167) 

involves a product of q components of the Q matrix and 

q-l q-1 

7e(N2B)' ;=0 '-'■■'^ j=Q 



is the [or] -dependent factor. In (II68I 1, the new notations 1 [fi] - fS and i [fi] - ft for 
a directed bond /3 € N2B have been introduced. These resuhs enable us to rewrite the 
mean field autocorrelation functions ( |164t as 

'^M = Z^^('^^ Z Fia](t,cr)IAt,cr), (169) 

where, for any t = (fo, . ■ • , f^-i) in (t, i)*, and for any cr e S q, 

I„(t,cr) = \im f fdQe-'^^7i(t,cr). (170) 

€^o4B(q- 1)! J 

The superintegrals /^(f, cr) are mean field superintegrals in the e ^ regime, and 
their values depend on the symmetry class. In the unitary symmetry class, the measure 

dQ reads 113 

1 dA\dA 

in terms of Efetov's polar coordinates, where Ai = cosh Q\ and A = cos 6, and the mean 
field source-free action S^^ is 

S^^ ^ 2B£{Ai - A). (172) 



dQ — ^^ , — —driidri\dKldK\d(pd^ (171) 



It can be checked, and it is stated in ||36l . that, in the unitary mean field superintegral 
( II7OI 1. the lowest order term in e is obtained by only retaining in ji(t, cr) its terms of 
highest order in Ai and by replacing the expression (Ai - A) in ( 11711 ) and illli with Ai. 
Therefore, in the expressions 

Qrr - UlBBUlB COSh0BU2BUlBB + U\BFU2F COS 0fU2FUlFB - '^ 

Qaa - VIBBV2B COSh 9bV2BV\BB + V\BFV2F COS 6fV2FV\FB - 1 

Gffl - -UiBBU2BSmh0BV2BVlBB + UlBFU2Fism0fV2FVlFB 

Qar - -VlBBV2BSinh0BU2BUlBB + VlBFV2Fisin0fU2FUlFB (173) 

for the components of Q^, which follow from ( 1155b . ( 1156b and ( 1157b . only the first 
terms of the right-hand sides contribute in the limit ( 1170b . Moreover, for the same 
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reason one can replace sinh 9i with cosh 6\ . These remarks, together with the formulae 
(fBST i, ( fT59l ) and ( fT60l l. lead to 



with 



and 



X s A] 






-\r]\\K\ 



U2B = 



RA 







Mcr, 



(174) 
(175) 



RA 



\r]\ = l-2r]l7]u \j<\ = 1+2kIki (176) 

In ( |174t and henceforth, for a and ii two functions of the Efetov polar coordinates, the 
equivalence a ~ b means that b can be substituted for a in the integrand of the mean 
field integral without modifying the result. The equivalence ( 1174b implies that 7T(t, cr) 
in (11671) satisfies 



7T(t, a) 



|^|2(,-l) 



v-?-! 



q-1 

11 '^''""l x„rx,„xt^ ifo-(0 = 0,/eN^-i 



'^aa-^rr 



ifcr(0) = 



q-\ 






?-l 



2^{q - 1) ■ KiK\ri\T]i ■ '^\W S,j,,^, 



(177) 



The last equivalence expresses the fact that only the term containing all the anticom- 
muting parameters can contribute to the superintegral ( 11701 ). Combining ( |170t . ( |171| l 
and ( 11721 ) with {Ai - /I) — > /li, and ( 11771 ) together, one arrives at 

4(<,cr) = lim e-^'^"^'A'f^-dAi-Yl6,.,,,^ n^/f,v(178) 

e^Q 2B(q - 2)1 Jo ^ l._l ''"'" {2B)i 1.1 ''"''' 

In the orthogonal symmetry class, a similar calculation leads to I„(t, cr) - {2B)^^. This 
result can also be inferred from ( 1169b . ( 11781 ), and from the expectation that the mean 
field intensity correlation matrix should satisfy C^^, - C^^, and should not depend on 
the symmetry class if a and a' are supported on two different bonds. 

By ( 11691 ) and by the results found above for I„{t, cr), the mean field autocorrelation 
functions become 






(179) 



(elT.i)' 



In order to get some explicit formulae out of (1179b . one can for example apply the rule 
( |109b once more, and notice that 



C 



[a] 



(2Byi 



(q-i 

n 

\k=0 



_d_ 



?-i 



det 



1 - Yj J'^^^"^'^ 



(180) 



*=o 



7=0 
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where, for j8 € N2B a directed bond, Nip) is the matrix acting on y[ i8> TR defined by 

N(J3)^\ lE^fi EP-P\ ^^^ (181) 



Let us first consider the unitary symmetry class ([/), and let us characterize the list 
[or] of directed bonds by another list j8 = {ji\,. . . ,/3„) of distinct directed bonds and a 
vector of integers q = {q\, . . .,q„) such that ySy occurs exactly ^y times in [or]. With the 
notation C^^{q) - Cj^f , an explicit calculation of the determinant in ( II8OI 1 shows that 



c-(.)-^n— ^ -fi 



.^0 ,.:W^ 



■' (182) 



Let us now consider the orthogonal symmetry class ((9), and let us characterize the list 
[or] by the list )8 of distinct directed bonds and the vector q of integers 

p = (^„...,A„,A„.:,...,A,) p.-(/^4) ,,,N„, (183) 

\ q^{q^,...,q„„q„,^u---,qn) [qj = iqj,qj) 

such that the components of q indicate the number of occurrences of the corresponding 
elements in y8. A first inspection of the formula (1180b shows that C^^{q) factorizes as 

m n 

Cfiq) = n ^ZS'lx^ qt) n Cffiq^) (184) 



k=l k=m+l 



Then, a calculation of the determinant shows that the correlation functions C¥^{qk) are 
given by the unitary formula (1182b . and 



-,MF 



c;i(^'§) = 



1 5« 5« 1 



{2B)'}+i dji d]i l-j-j 



j^ho W'^" 



(q + qy. 

^^ ^' (185) 



It can be checked that these formulae coincide precisely with the predictions of the 
Gaussian Random Waves Models (08]) and ( |53] |. 



Moreover, it can be checked in (11791) that the autocorrelation functions C7, defined 

at the end of Subsection l4. 1 1 give rise to mean field autocorrelation functions C'^\ that 
do not depend on the particular convention cr e S q chosen. 



6 The Gaussian Correction 
6.1 Beyond Mean Field Theory 

It is known IfTTI that not all increasing sequences of quantum graphs are quantum 
ergodic, and hence, the mean field theory does not always yield the main contributions 
to the autocorrelation functions. Therefore, it is necessary to estimate the importance 
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of the supermatrices Z and Z lying off the zero mode manifold. For this purpose, let us 
write 

Z B Zo + 6Z, Z = Za + 5Z, (186) 

with (Zo,Zo) on the zero mode manifold and with {5Z,5Z) an orthogonal deviation, 
and let us expand the exact action S [Z, Z] in (|149t up to second order in 6Z and 6Z 
around (Zq, Zq). The truncated action S [SZ, SZ, Zq, Zq] obtained in this way leads to a 
generating function 

diZo, Zo) d{6Z, 6Z) e'^ I-^Z-^z^Zo.Zo] ^ ^ j gy^ 



/' 



Suppose for the moment that Z and Z in (1187b are also required to satisfy Z = Z^. It 
follows that the partial traces over ^ of the supermatrices 6Z and 6Z in ( |186t must 
vanish. This property implies that the truncated action S has no linear terms in 6Z and 
6Z. A direct but tedious calculation shows that, if the sources in the truncated action 
are set to zero, one gets 

S [6Z, 6Z; Zo, Zq] = 5 ^ fZo, Zq] + S <"' [W; w] , (188) 

where S^^[Zq,Zo] is the mean field source-free action, 

W = ( 1 - ZoZo)" ^-6Z{l-ZoZo)-K W s ( 1 - ZoZo)" MZ( 1 - ZoZq)" K ( 1 89) 
and Sq [W, W] is the term of the exact source-free action S'oIW, W] of second order in 



W and W around the origin. If e^^" was to be integrated as in (11871 ), the changes of vari- 
ables ( |189t , which both have unit Jacobian, would factorize the superintegrals over the 
zero mode (Zq, Zq) and over the orthogonal deviation (W, W). This factorization occurs 
because the domain of the superintegral over {W, W) is independent of (Zo, Zo). Indeed, 
it can be readily seen that the equality Z - Z' merely becomes W = W^. Moreover, 
if all the Grassmann generators of A are sent to zero, which is what really matters for 
the domain of a superintegral, the color-flavor requirements ( I144l i on (6Z, 6Z) become 
Wbb - Wgg and Wff - -W'j,^, and there is no further condition concerning the eigen- 
values of the positive Hermitian matrix W^gWBB- Finally, the condition ensuring that 
(6Z, 6Z) is orthogonal to the zero mode manifold forces W and W to have vanishing 
partial traces over ^. 

By analogy with the situation described above, where the sources are set to zero 
and where the supermatrix variables are constrained to satisfy Z - Z^, one defines the 
truncated generating functions 

lwC/) = CC/)-f&]C/), (190) 

where the Gaussian generating function is defined by 

d''(W,W)e-'"'^'^-'^^ (191) 



€)0')^J^ 



and S^^^[W, W] is the term of the exact action S [W, W] in ( |149t of second order in W 
and W around the origin, namely 

S^-\W, W] = str| WW - ^J',W^J'«W - -WS'lWSA . (192) 
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The integration in ( |191| l is over all supermatrices W and W in L{TR ® J[\TR ® ^) 
that are diagonal in ^h, define a configuration (W, W) orthogonal to the zero mode 
manifold, and satisfy the color-flavor conditions Wbb - W^^ and Wfp - -Wp^. The 
measure d'^{W, W) is then the product of the flat Berezin measures over the independent 
components of W and W. 

In ( 1112b . it is claimed that if either the advanced or the retarded sources are sent to 
zero, the exact generating function becomes identically equal to one in a neighborhood 
of the origin. It can be checked that the same property holds separately for the mean 
field and the Gaussian generating functions. Therefore, if the formula dllOt is used to 
define truncated autocorrelation functions €[„] from ^[a], one gets 

Cw=C^]+C[;,, (193) 

where C^ are the mean field autocorrelation functions found in the previous section, 
and 



(26)''-' ... ... (f\ d^ 



^'«i=l"S2^(^rT)7<«i' '€^^ 



n.. 



^m(O), (194) 



are the Gaussian autocorrelation functions. 

In fact, in order to calculate the Gaussian generating functions ( 11911 ), one can first 
calculate the second order generating function 

^'^^U)^jd'^'\Z,Z)e-''"^^-^\ (195) 

defined from (I191l l by relaxing the constraint that (Z, Z) must be orthogonal to the zero 
mode, and then divide by the second order mean field generating function 

C''^U)^ f d^^^'\Y,Y) e-'"''-'^''^\ (196) 

which contains the zero mode contribution to ( 1195b . With the notations Zq = Ij^ ® 7 
and Zq - lj?i ® Y for the supermatrix variables in the zero mode, d^^^^\Y, Y) is the 
measure induced by S^\Z,Z) on the zero mode manifold. Similarly, the mean field 
second order action S^^^~^ is obtained by restricting 5*"', that is 



5MF(2)[j, y^ ^ 1 jj.^(2 - e-2^)l^ ®YY- JrYjj) . 



(197) 



The generating functions ( |195b and ( |196b are identically equal to one in a neighborhood 
of the origin if either the advanced or the retarded sources are set to zero. Therefore, 
the Gaussian autocorrelation functions can be written 

with the obvious definitions for SC \ and S^, / . 
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6.2 Diagonal Modes in Direction Space 

Let us first only consider the subset of supermatrix variables (Z, Z) that are diagonal in 
the whole amplitude space J{ - Jii, ®J{d. If d, J JZ,Z) denotes the measure induced 
from (P'\Z,Z) on this subset, the goal is to calculate the {dd) second order generating 
functions 

Cl(..)C/-)-/<^:.)(2'Z).--t-.^l (199) 

The diagonal modes Z and Z are parametrized in time-reversal space as follows, 

and their generalized transposes read 

The modes Z, , ZJ , Z3 and Z3 are only considered if time-reversal invariance is con- 
served, so that Z, Z and their generalized transposes all become diagonal in time- 
reversal space. When these formulae are substituted into 5*^'[Z, Z] given by ( |192| i. 
the diagonal modes in time-reversal space, which are indexed by 1 and 4, are coupled 
together, and do not mix with the off-diagonal ones indexed by 2 and 3. After some 
algebra, one finds 

S%,-Sf]^^+{K-l)S'^^^ (202) 

with the diffusion action S 'J^ and the cooperon action S .L defined by 

SfJI^^ = str(ZiZi+z]zJ-y,z]7,Z, -Zi5^z|5,) (203) 

Sf]^^ = str(z|z]"+Z3Z3- 7,237^2! -Zl^r'Za^,) (204) 

Hence, the generating functions C„i i^^-, factorize as 

A2) ^A2)D (A2)C y-i ,205) 



where, for o e {D, C), 



€udd,U)-ld^iZ,Z)e-^^^''K (206) 



In ( I2O6I 1, the diffusion measure is the product of the Berezin measures of the inde- 
pendent components in the supermatrix variables Zi, Zj, Zi and z], and the cooperon 
measure is similarly formed with the independent components in Z,, Z3, Zj and Z3. 

Notice that the cooperon generating functions (I206l l, which only exist if 5^ = S , 
can be obtained from the corresponding diffusion generating functions by replacing 
Ja with J^ . One can thus focus on the diffusion generating functions and infer the 
cooperon result from this remark. 
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In order to perform the diffusion and cooperon integrals in (12061 ), the supertraces 
in ( 1203b and ( 1204b must be explicitly expanded in Bose-Fermi space. Notice that the 
color-flavor conditions on the Bose-Bose and Fermi-Fermi components of Z and Z in 

( |200b become Zjbb - ^bb ~ ^*bb ^^^ ^Jff ~ ~^ ff ~ '^*ff ^^'^ ^ J ^ ^4- ^^^ ^^ 
define for each directed bond /3 e Nib 



,t _ 



--{z; 



IBB' ^\FF' 



ABB' 



"4FF 



)P 



-\QI3 



1^ 
40/3 




(207) 



The vectors ZiQ/j ^"d Ziq^ contain the commuting parameters of Zi^ and Z\p respectively, 
and the vectors z - and z - contain those of zL and zL. Similarly, the anticommuting 
variables of the diffusion action are arranged in the vectors 



ini3 - (^iBF , Zifb) 



^4l/3 ~ ^^4BF ' ^AFB^ 



41/3 



ZlBF 
Z\FB 

_4BF 

ryif. 

^AFB 



(208) 



Collecting the IB row-vectors Z\^p (resp. z^, ) together, one can write a larger row 
vector Zio (resp. z' ,). The column-vectors Z\^ and zl= are formed similarly from Zj^g 
and z' - , and one proceeds in the same way with the anticommuting variables in ( |208b . 
Let us also introduce a 2B x2B matrix s defined from the Bose-Bose blocks of the 
source supermatrices Ja and J, by 



s(j)m = JaUa)BB,P'l3J rUr)BB,pp' ■ 

A direct expansion of ( 12031 ) in Bose-Fermi space then leads to 



c{2)D _ .= _ J J 1:/I®l2x2 -( 



A2)D 



s(0,0) 
lj?l ® l2x2 



for the part o^ S ~L involving the commuting variables, and 



,(2)D 
¥d)l 



(2il > z\^ 



M, ®l2x2 lj^®(-ii) 



^10 

z"'"- 

^40 



Zll 

2!\ 



(209) 



(210) 



(211) 



for the part involving the anticommuting variables. Notice that these formulae depend 
on the scattering matrix S^- e^^S only through the classical map M^ = e^^'^M it 
generates. It is straightforward to calculate the diffusion superintegral ( 12061 ) from these 
quadratic forms, and the cooperon generating functions are found by substituting the 
matrix f^ for the matrix Ja in the diffusion results. If o e [D, C), one gets 



(2)„ . . det(l^-,9°(;„,0)M,)det(l^-i°(0,jV)M,) 



det(l,^ - M,)det(l^ - s°(jaJr)M,) 



(212) 
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where s^ s s in ( 12091 ). and s*- is obtained from s by replacing Jg with 7,f . 

In order to unveil the mean field contribution to (12121 ). one restricts the superintegral 
in (|199t to the zero mode Zq = lj?i ® Y, Zq = 1^ ® ?. The supermatrices Y and Y are 
required to satisfy Y - Y'^, and can thus be parametrized by 

Yo Yc\ ^ y I Yu <y^ 



Y-\ vr^BF vf and y = U V^ (213) 



in time-reversal space. Then, the second order mean field action ( |197t splits into a 
diffusion part, containing the supermatrices I'd and Fp, and a cooperon part, involving 
Fc and ?c- One gets 

5MF(2)D ^ str((2-e-2.^FD?D®l.^-/,i'D^?D) (214) 

5MF(2)c ^ ^xi[[2-e-^')YcYc®ty^-JrYcJlYc). (215) 

It follows that the mean field contribution to (1212b factorizes into a diffusion and a 
cooperon factor, as in ( 1205b . The expressions (1214b can be developed in Bose-Fermi 
space, and the resulting quadratic forms have inverse superdeterminant 



,MF(2K. .^_ (l-^-^--^°(^0))(l-^-^--^°(0, jV)) 



C^ U-^Jr)-- -, — -^ '-, °e{Z),C), (216) 



where 



(217) 



0-°(ja,jr) = J^ Yj ■^°(./«'^''-W ~ 1- 

Notice that no index (dd) has been added in the left-hand side of ( 1216b . The reason 
is that the zero mode supermatrices Zq and Zq are always diagonal in the direction 
space y[^, and hence ( 12161 ) is also the mean field contribution to the full second order 
generating functions ( 11951 ). 

The diffusion and cooperon generating functions in ( 1212b and ( 12161 ) become identi- 
cally one in a neighborhood of the origin if either the advanced or the retarded sources 
are set to zero. Hence, the (dd) Gaussian autocorrelation functions defined as in ( |198b 
from the (dd) second order generating functions ( 1205b split 

^[alidd) - ^[aUdd) + ('^ ~ l)^[«],Wd)' (^18) 

where the diffusion and cooperon autocorrelation functions in the right-hand side are 
defined by the formula ( 11981 ) applied to the diffusion and cooperon versions of ( 1212b . 
that is 

^S,(.0 - i"S ^15^ («(..) - <n ' - 1^' C). (219) 

The next step towards the calculation of ( 12191 ) is to calculate the derivatives of 
the diffusion and cooperon (dd) second order generating functions ( 1212b . Performing 
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explicitly the unique advanced derivative, one easily finds that for all integers q > 2, 
and for all sets of q directed bonds [or]. 



se 



.(2)o 
aUdd) 






1 



|:(«'^'^>^T-Ao,,v)M. 



(220) 



7V=0 



where jV = iji, ■ ■ ■, iq-\Y are the q- I retarded sources. 

Let us first consider the situation q -2 and [or] = {a, a']. Then, the unique retarded 
derivative in (12201) can be performed, and one gets 



5^^'^° = tr 



M, 



M, 



M, 



s„ + 



1 - M ' 1 - M '' 1 - M, 



(221) 



where s°, 5° and 5°^ respectively denote the derivatives of 5° with respect to jr, ja, and 
jr and ja, all evaluated at yV = ja - 0. If these derivatives are calculated using the 
parallel convention for Jg and Jr, the diffusion (dd) derivatives take the form 



J2)D 






M, 



1 -Ms 



M, 



M, 



1-MJ,„,\1-M,,„,, 



(222) 



and the cooperon (dd) derivatives is given by the same formula with a in place of a. 
Notice that this expression agrees with the derivatives ( I124l i and ( 11251 ) of the diagonal 
approximation to the trace formula for the generating function in Section 14.21 If the 
derivatives of s° found with the crossed convention are plugged into ( 12211 ). one gets the 
diffusion (dd) derivatives 



"^{a,a'\ldd) 



M, 



1 -M, 



+ 5„ 



M, 



1 -M, 



M, 



aa \ ^ I eta 



(223) 



and the corresponding expression with a in place of a for the cooperon (dd) derivatives. 
This is again the result obtained in (1124b and ( 11251 ) for the derivatives of the diagonal 
approximation to the generating function. The derivatives of the second order mean 
field generating function ( 1216b have to be removed from the previous formulae. They 
read 

1 \' I \ 



c/;MF(2)o _ 



1 



g-2e 



0-,0-a + 



1 



0"ra, 



(224) 



where the indices r and a denote the derivatives taken on cr, which are all evaluated at 
the origin. These derivatives can be calculated according to the parallel or the crossed 
conventions. The results obtained are given by the formulae ( 1222b and ( 12231 ) by system- 
atically replacing the sum of classical walks Me( 1 - M^)" ' with the uniform component 
of (1 - Me)"', which is defined as in ( |35] | and reads (1 - e"^'^)"'|l)(l|. This draws a 
parallel between the zero mode and the uniform component. Finally, when the formula 
( 12 19b is applied, the terms in ( 12221 ) and ( 1223b that are too singular in e are exactly 
compensated by the mean field derivatives, and one is left with the finite result 
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a,a',{dd) 
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xG.D 

aa' Xdd) 
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QiY'Xdd) 



(225) 
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for the parallel and crossed diffusion {dd) Gaussian intensity correlation matrix. In 
these formulae, the matrix R denotes the massive component defined by the decom- 
position ( [35T l. The cooperon contributions to the Gaussian intensity correlation matrix 
read 

-,G,C _ R &a' + Rg'S, _ 2 ^^^ ^xc,C _ ^ r-G,C 

(2B)2 (2Bf 



^aa',(dd)~ (1R\1 (7 R^i ^ ^ aa' ,(dd) ~ ""•'^'^ aa' ,(dd) (226) 



The surprising second terms of the formulae in ( 1225b and (12261 ) originate from the 
fact that the zero mode contribution to (1222b and (12231 ) is the uniform component of 
(1 - Mf)"' and not M^{\ - M^Y^ . This discrepancy is due to the additional symmetry 
Y — Y'^ of the zero mode. Notice however that these second terms are of higher order 
in (2B) ', and are thus of minor importance when the large graph limit is considered. 

Using the same strategy as above, it is not difficult to calculate more retarded deriva- 
tives in ( |220b and to remove the mean field contributions. If [or] is a list of q' > 2 
directed bonds, and if the convention cr = id e 5^^ for the generating function ( 1212b is 
chosen, one gets the formula 

1 "'^ 

for the diffusion {dd) Gaussian autocorrelation function of degree q, and the same 
formula with ao in place of ffo for the cooperon {dd) Gaussian autocorrelation function 
of degree q. 

6.3 OfF-Diagonal Modes in Direction Space 

Let us now investigate the full second order generating functions taking into account the 
modes Z and Z that are off-diagonal in direction space. The parametrizations ( 1200b of Z 
and Z in time-reversal space can be kept, and hence the formulae ( 1203b and (1204b also 
hold in the presence of off-diagonal modes. This implies in particular that the second 
order generating functions factorize into diffusion and cooperon generating functions 
as in (1205b . and that the cooperon formulae, which are considered only if time-reversal 
invariance is conserved, can be found from their diffusion counterparts by replacing 
Ja{ja) with Ja{ja)''' ■ One can thus temporarily concentrate on the diffusion modes only. 
One can distinguish between diagonal and off-diagonal modes and introduce the 
notations 

Z''^''^ = Zpp and Zf = Zpi„ (228) 

and similarly for Z. The quadratic action couples diagonal modes with themselves, 
which is precisely the part treated in the previous subsection, off-diagonal modes with 
themselves, and diagonal modes with off-diagonal modes. 

The integration scheme used here is similar to the one that leads to the explicit 
formula (1212b for the {dd) Gaussian generating functions in terms of four determinants. 
Let us first focus on the commuting components Z^,,, s € {B, F], of the fields. The row 
and column vectors defined in ( 1207b . whose purpose is to write the diagonal action 
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S ,L as a quadratic form, are adapted to the situation where the fields Z and Z are as 
in (12281). Let us define 



/T. - /^diag tdiag -off _toff\ 
^/^-I^IO '-40 '%'^- ' 



40 Jb 



(229) 



where 2 "f^ and z,-!^^ are formed with the diagonal modes of Zia and z , defined in 
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( |207t . and 2°? and z'? are formed with the off-diagonal ones. We proceed in the same 
way with the column vectors and introduce 



W/j 



_ / diag T -tdiag T ^off T -toff r'i 
\~10 ' ~40/3 ' ~10 ' ~46/J j 



(230) 



Then, a careful inspection of the diffusion second order action ( |203l l and some algebra 
show that the part of this action involving the commuting variables is the quadratic 
form 
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defined by the 16B x 16B matrix 
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2x2 



l^i® 1 



(232) 



2x2 / 



Notice that the block coupling the diagonal modes together is precisely (I210l i. In the 
quadratic form ( 1232b . the 2B x 2B matrices P^, Q^ and K^ are defined by 



Pepp' = S epp'S 



epp" 



Qepp' = Sepp'S*. and K^/^/s, = feo,5^„s5 



epp 



P.P'^ epp^ efSp- 



(233) 



In fact, Pj and Q^ both vanish since we only consider simple graphs. The square of 
K = lim^^o K^ is the diagonal matrix 



{K\^5pp,MppM-pp, 



(234) 



which only depends on S through the classical map M. It can be deduced from (I234l l 
that the spectrum of K is real and is contained in (-1, 1) if the graph is ergodic. In 
(1232b . s{j) is the matrix defined in (1209b . c{j) is another matrix satisfying c(0) - 1, and 
a{j) and b{j) are given by 



and 



aUa,jr)pp' = Spp'jaE^I^I^ + ^0, [jrE^'^)^ 
b(jajr)pp' ^ SppJ^E'^'l + dpp. {jrE^"-') 



(235) 



(236) 



pp 'PP V Jpp 

It can be checked that a(j) and b{j) both vanish if the convention o" = id for the 
generating functions is used. 
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The determinant of S in ( 12321 ) can be calculated using the rule ( I134l l adapted to 
conventional determinants. The result is the product 

detS = Aff-ff ■ Abb-bbU) (237) 

where Aff-ff = Abb-bb{0), 

Abb-bbU) = det (1:^ - sU)M,) det (1^ - cUW,) det (1^ - N,U)) , (238) 

and 

1 - sC/)A^£ 1 - c(j)K, 



The first factor in (12371 ) comes from the couplings between FF and FF components 
of the variables in the vectors vv/j and w^. The second factor comes from the couplings 
between BB and BB components. The fact that the contribution of the FF - FF cou- 
plings can be found from the contribution of the BB - BB couplings by setting all the 
sources to zero can actually already be observed on the formula ( 1203b for the diffusion 
second order action S '^*^. A further look at this formula enables one to deduce the con- 
tributions of the couplings between the anticommuting variables. It can be seen that 
the matrix mixing the BF components of the row vectors vv^ and the FB components 
of column vectors w^ has determinant ABB-BBija,^), and similarly, the matrix mixing 
the FB components of the row vectors vv^ and the BF components of the column vec- 
tors Wj3 has determinant ABB-BB{0,jr)- Hence, the diffusion second order generating 
function reads 

A2)D ^ A2)D A2)D ADD ,r..r.. 

'^[a] '^[a\,(dd) ^[o-],(oo) '^[a\,(doy \^^^J 

where the first factor in the right-hand side is the diffusion {dd) second order generating 
function (1212b . 

(2)^ . . _ det (l^ - c{ja, Q)K,) det (l^ - c{QJ,)K,) 
' '-^ ^ det (l^ - K,) det (l^ - c{ja,jr)K,) 

and 

^^^^ . . _ det (l^ - N,Ua, 0)) det (l^ - N.iOJr)) 

det [t^ - N,{Q, 0)) det (l.^ - N.ijajr)) 

These functions all have the property that they become identically one in a neigh- 
borhood of the origin if either the advanced or the retarded derivatives are set to zero. 
Hence, their product ( |240b . and the cooperon analogs, share the same property. It fol- 
lows that the derivatives of these functions satisfy 

<- Z ^^S!v) + (--l) Z <iur (243) 

x£{dd,oo,do\ x£{dd,oo,do} 

Moreover, since K has no eigenvalue unity, the (do) generating functions at e = are 
analytic in a neighborhood of the origin, and hence, their derivatives cannot contribute 
to the Gaussian autocorrelation functions ( 1198b . 
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If the convention cr = id is used for the generating functions, then the functions 
a and b in ( I235l l and ( 1236b vanish. In this case, the (do) generating function ( I242l l 
is equal to one, its derivatives vanish, and only the (dd) derivatives remain in (|243l l. 
Therefore, from ( |243t , (|227t and this remark, the Gaussian autocorrelation functions 
( |198l l of degree q read 



*#/ 

Here, the last term {cq — > Q'o) stands for the right-hand side with ao replaced with 
ao- In summary, the ofF-diagonal modes do not bring any additional contribution to the 
Gaussian autocorrelation functions if the convention cr = id is used. 

Let us now consider the situation where q - 2, [or] = [cc, a'], and the generating 
function is defined with the crossed convention. In this case, the two derivatives on the 
{do) generating function ( 1242b give 



<Jff^' n M . = tr 



1 1 1 



1 - N, • l-N, 1 - A^, 



(245) 



where A^^ denotes the value of the function N^ijo, j\) at the origin, and A^^ o, A^^ i and 
N^fi] stand for its derivatives at the origin. Since a and b are zero at the origin, A^^, N^fi 
and Nej vanish. Therefore, only the second term in the trace of ( 1245b contributes, and a 
short calculation shows that the diffusion (do) Gaussian approximation to the intensity 
correlation matrix in the crossed convention reads 



aa' ,(do) ^^0 g "-"Xl-M, /,„. V\\-Klaa \1-KI. 



(246) 



The cooperon result turns out to be the same. With the decomposition ( |35] | of the 
classical walks, and using the fact that the diagonal elements of K" vanish if the integer 
n is odd, one gets 

where 

The matrix R'^, called the back-scattering matrix, is formed with all the oriented walks 
followed with the classical map M which involve only back-scatterings/? -^ p and no 
transmission. Together with ( |225b and (|226b . ( 12471 ) yields the Gaussian contribution 

C^tt' = 7Tm2 p^ff» + {.K- l){Rda + Raa)\ - ^a.a' 



(2B)2L "" ^ '^ "" ""^J "•" {2BY 



-\-K- 



-(2fi)2K+^y (249) 

to the intensity correlation matrix in the crossed convention. 

For q - 2 and a - a', the parallel and crossed results ( 12441 ) and ( |249b coincide. 
Notice that in this case, the parallel and crossed sums over oriented walks represented 
in Figure|2]are also the same. 
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7 Criteria and Rates of Universality 

7.1 Full Universality and Criterion for Ergodicity 

The calculation scheme summarized in ( |193t leads to truncated autocorrelation func- 
tions, which are sums of the mean field contributions obtained in Section |5] and the 
Gaussian contributions obtained in Section |6] On one hand, the mean field results co- 
incide with the universal Gaussian Random Wave Models introduced in Subsection |3.2| 
On the other hand, the Gaussian quantities in (1244b depend on the quantum graph, but 
do so only through its classical map M, and more precisely through its matrix R, which 
is defined in (35) and represents the massive component of the sum of classical paths 
■j^-. The importance of the truncated autocorrelation functions is twofold. Firstly, 
their Gaussian contributions can be compared with their universal mean field parts in 
the limit of large graphs. These comparisons lead to conditions on the increasing se- 
quence of quantum graphs {(G/, 5'/));eN to asymptotically follow the predictions of the 
universal Gaussian Random Wave Models, in which case we say that full universality 
is met. Secondly, for a class of increasing sequences larger than this universal class, 
the truncated quantities approximate the exact autocorrelation functions (l42l i. 

Let {G,S) be an ergodic simple quantum graph, and let [a] be a list of q directed 
bonds for some integer q > 2. The mean field autocorrelation functions found in 
Subsection l5.2l and the Gaussian autocorrelation functions in ( I244l i lead to the truncated 
autocorrelation functions 

^ c_OfAa]^ ^ 1 ^ r ^ I ^ ji _ _J^ 5 

yi 

The parameter k is equal to one or two depending whether time-reversal symmetry is 
broken or conserved, and c{k, [a\) is the combinatorial factor defined in the following 
way. Suppose that each directed bond/3 appears exactly qp times in the list [or], then 
c(l, [a]) = Hflf 1 1/3- Suppose now that pi, denotes the number of directed bonds in 
the list [a] supported on bond b, then c(2, [a]) = Y\b=i Pb- The intensity correlation 
matrix and the moments implied by ( I250I I read 

and 



1 K- 1 

—J (1 + Saa' + Raa' + Ra'a) + . p2 ^^aa' "*" ^aa' "^ R&'a) (251) 



^?.'>- = 7^ 1 1 + 7 TTT^'" + '^('^"Q + t^&a + iq- 2)Raa)] ■ (252) 

ilBY'X (q-iy. q\ I 

In these two expressions, the last term of ( |250t , which is of higher order in the inverse 
number of bonds B"', has been neglected since we are ultimately interested in the large 
graph limit B — » oo. 

For graphs in the orthogonal class, the exact autocorrelation functions C\a\ in ( l42l l 
do reflect the symmetry |fl„|- = la^^p of the wave function intensities. The intensity 
correlation matrix ( 125 11 1 indeed satisfies Caa' - Ccm' if k - 2, but the truncated auto- 
correlation functions ( I252l i of degree ^ > 3 do not obey such a symmetry in general. A 
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comparison of ( 125 11 1 with a numerically obtained intensity correlation matrix for com- 
plete quantum graphs with various choices of scattering matrices reveals that ( 1251b 
captures the asymptotics B ^ oo very well at least as long as the intensity matrix is 
not dominated by the massive contributions (see Figure[3j. A numerical comparison of 
higher moments to (|252t (not shown) reveals that their massive contributions are not 
as well approximated by our theory. We should however emphasize that the numerical 
evaluation of higher moments is not very stable. The statements we put forward below 
are all consistent with the numerically obtained data we have. 

Let {(G/, 5;))/eN be a sequence of increasing ergodic simple graphs that are either 
all in the unitary or in the orthogonal class. The formula ( |250l l suggests that the mean 
field term is the dominant one if and only if the sequence of matrices {^/)/eN converges 
to zero as / — » oo. Consequently, we conjecture that the Gaussian Random Wave 
Models ( |48] | or ( l53b are asymptotically met in the sequence {(G/,5'/))/eN if and only 
if the norm of the matrices Ri decay as Z — > oo. Note that Ri —^ Q introduces a small 
parameter in which a systematic expansion may be performed. We strongly believe that 
our approach may be extended to a rigorous proof that Ri ^ Q implies convergence to 
the Gaussian Random Wave Model. 

Notice that, by the definition ( 1351 ) of the matrix R^, the decay of the sequence R^j 
is equivalent to 

lim — ^ = 0(6). (253) 

This is also equivalent to the decay of Ri = lim^^o R^j. Since all the components of the 
matrix M^j = e^^'^Mi are non-negative, ( I253I I implies that 

lim Ml = 0. (254) 

A necessary condition for this property to occur, and hence for full universality to be 
met, is that the valencies of the vertices all tend to infinity. This condition has already 
been derived at the end of Subsection 13. 21 where the Random Wave Models are built. 
Conversely, suppose that (1254b is fulfilled, then the equality (1253b also holds, and thus 
Ri decays. Therefore, the necessary and sufficient condition for full universality that is 
conjectured above is actually equivalent to ( 1254b . 

The expression (1251b . together with the formula ( |69] l, generate a truncated expres- 
sion T'v for the fluctuations of an observable V. The asymptotic quantum ergodicity 
problem described in Subsection l3.3l can be addressed in terms of these truncated fluc- 
tuations. Moreover, in the situations where 'Fy decays, 'Fv is expected to approximate 
the exact fluctuations Fv- A direct calculation shows that, for an observable V with 

^ tr(yL)^ ,, i:,,'[VL-R-VL]^,, 

This formula motivates the following criterion for asymptotic quantum ergodicity to 
be met in an increasing sequence of quantum graphs. We conjecture that an increasing 
sequence {(G/, 5/))/gN of ergodic simple graphs is asymptotically quantum ergodic if 
and only if 

}'^ (^^ ^ ° ^^''^ 

51 







Figure 3: Rescaled matrix elements of the intensity correlation matrix for complete 
graphs with V vertices. 

The four panels correspond to four different choices of the scattering matrix. The upper 
left panel is a complete DFT graph. For the upper right and lower left panel the scat- 
tering matrix that corresponds to one vertex has been chosen by the following unitary 
transformation of a DFT matrix o-'{A) = e^^^'^^^^^^S^rr^rw with the values A = 0.5 
for the upper right panel and A - 0.25 in the lower left panel. For the lower right panel 
Neumann scattering matrices have been chosen. 

An additional magnetic field was applied to break time-reversal symmetry. Black sym- 
bols and lines correspond to results for the orthogonal class and blue symbols and lines 
to results for the unitary class where the symbols correspond to numerically obtained 
intensity correlations and full lines to the corresponding prediction ( |251t . The squares 
(orthogonal) and downwards pointing triangles (unitary class) give the average rescaled 
diagonal element 2B Yja=i ^»« ^^^^ Gaussian Random Wave Model predicts the value 
2 indicated by the upper dashed line). The diamonds (orthogonal) and upwards point- 
ing triangles (unitary class) give the average rescaled time-reversed diagonal element 
2B Yj7it\ Caa (the Gaussian Random Wave Model predicts the values 1 in the unitary 
and 2 in the orthogonal class). Note that the diamonds and squares always lie on top 
of each other Eventually, crosses (orthogonal) and circles (unitary class) give the av- 
erage rescaled element -^ TH^ a=i-ai=a' a^a' ^aa' ('^^e Gaussian Random Wave Model 
predicts the value 1). The corresponding predictions from ( 1251b are given by the full 
lines. 
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for any acceptable sequence {V/I/en with V; - 0. We will give a slightly more detailed 
variant of this conjecture (and a discussion of possible obstruction to its validity) below. 

Moreover, if the stronger condition /?/ ^ is fulfilled, the increasing sequence 
of graphs is fully universal, and the convergence rate of 'Fv, and hence of 'Fv, is then 
also universal. As in the case of the validity of the Gaussian Random Wave Model we 
strongly believe that our approach can be extended to a rigorous proof using R (or an 
equivalent quantity) as the small parameter. 

Note that the crossed formulae for the massive (Gaussian) contribution to the in- 
tensity correlation matrix in ( |249t differ from the formulae in ( 12441 ). which are used 
above. The crossed expressions only involve the diagonal components {a, a) and the 
components {a, a) of the matrix R - they also contain a new backscattering term. These 
do not obey Caa - Caa and, indeed, they do not capture the massive corrections in the 
exact correlation matrix as well as the parallel convention in a numerical test (Figure[3] 
only presents the results for the parallel convention). 

Let us now consider observables V such that VtLt, = § on half of the bonds, and 
VbLf, = ~|§ on the other half. The set of such observables is actually sufficiently large 
to compare the intensities of the wave function on the different bonds. Moreover, they 
provide acceptable sequences, according to ( |63] |. For such observables, ( 1255b yields 

1 2Ktx{R + Rcri) 

Tv ~ K— + ;-^, (257) 

2B (2B)2 ^ ' 

where we neglected almost all off-diagonal terms (apart from those obeying a' - a) in 
the double sum in ( I255I I. An increasing sequence {{Gi,Si)]im of simple graphs is then 
expected to be asymptotically quantum ergodic if and only if both 

,^llL^^^^f\zJ!!l±^Q (258) 

M« (2B,)2 /^oo {2Bif j^ mii 

and 

tr/^/crf 

hold. In ( 1258b . the complex numbers m/ ,, 2 < / < 2B, are the 2B - I non-zero masses, 
that is the 28 - I eigenvalues of the matrix 1 - Mi. 

7.2 Quantum Ergodicity and the Classical Spectral Gap 

Sufficient conditions for the condition ( 1258b to be fulfilled or violated in an increasing 
sequence {{Gi,S i)}ien of ergodic simple graphs can be given in terms of the sequence 
{AM,}zeN of spectral gaps of 1 - M/. 

Let us first consider the case that all non-zero masses stay away from the origin. 
The sum in ( 1258b behaves like 2B;, and hence, after dividing by {2Bi)^, the large graphs 
limit vanishes, and ( 1258b (and similarly ( |259b ) holds. 

Now let us turn to the case that some masses approach zero as Bi — » oo. For sake 
of simplicity, the index / of the quantum graph {Gi,S i) will be dropped. Let us order 
the spectrum {OT,),EN,g of 1 - M such that |m/| < |m,+i| for all / € N2B, and let us now 
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suppose that the spectral gap A^ = |'«2| approaches the origin with an exponential rate 
a > 0, that is 

Nzl ~ {2B)-'\ (260) 

The matrix R is real since M is real and the vector |1) is also real. It follows that the 
massive contribution of the fluctuations ( I257l i can be written 

^M _ 2/c sr^ ^-mi _ 2k ^ "Ktni 2k{2B - 1) 

^y - (2B)2 Z^ "^ mi (2B)2 Z^ |m,p (2B)2 ' ^^'''^ 



The second term of the right-hand side behaves like (2B) ', so that (1258b is satisfied if 

V 



and only if the first term of the right-hand side, denoted by 'F^ in what follows, decays. 



With the obvious inequality 5lm, < |m,|, one gets 



2k ^ \ 2k 2B-\ „ , 

(2B)^ f^ |ot;| (2B)- |m2| 



Therefore, if a < 1, 'F^ decays and ( I258l l is fulfilled. 

Since there are 2B - I non-zero masses, and since these masses are either real or 
appear in complex conjugated pairs, there is at least one mass m/ , such that 

ti = sup tanargm/,, < oo. (263) 

Remember that the non-zero masses are confined in the open disc of radius 1 centered 
at 1. It follows that |ot,| < (1 + f,)*Xm;, and thus 



-M ^ 2 %mi ^ 2 1 ,Tcvv-2 

^ (2B)2 |m,|2 - (2B)2 (1 + f,)N,| 



^f > TTT^ ^ ^ T^TT^ ., . .., , ~ (2B)"-^ (264) 



Therefore, if a > 2, F^ does not decay, ( 1258b is violated, and the increasing sequence 
of quantum graphs is not asymptotically quantum ergodic. 

7.3 The Four Possible Regimes 

To summarise our findings and give a more detailed account of our conjecture and 
possible obstructions to its validity we have found the following four regimes 

1 . The fully universal regime Ri — > 0. Equivalently, all matrix elements of the clas- 
sical map converge to zero, or all eigenvalues apart from the (non-degenerate) 
eigenvalue one of the classical map Mi converge to zero, or all 2B - I non-zero 
masses converge to one. In this case the eigenvalues of the classical map may be 
used as a small parameter for a systematic expansion. Our theory then shows that 
all autocorrelation functions converge to the universal predictions of the Gaus- 
sian Random Wave Model, and we believe that the scheme used here can be 
extended to a rigorous proof. Such graphs are also asymptotically quantum er- 
godic with a universal rate of convergence. 

This class includes the complete DFT graphs (or complete quantum graphs such 
that nonvanishing elements of the classical map are of order C)(y ')). 
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2. The large mass regime, characterized by a < 1 in ( 1260b and Hm/^oo Ri + 0. Some 
non-zero masses do not converge to one, but none of them approaches zero too 
fast. The complete Neumann graph is an example in this class with a - \. Our 
theory predicts massive corrections to the predictions of the Gaussian Random 
Wave Model which persist in the asymptotic regime B — > oo. This implies that 
the Gaussian Random Wave Model is not applicable but asymptotic quantum 
ergodicity still holds. This statement has the status of a conjecture which may 
be very hard to prove because there is no obvious small parameter As a con- 
sequence our result ( 1250b may not estimate the massive contribution accurately. 
Our numerical data (see Figure[3]l indeed show that the intensity correlation ma- 
trix for a complete Neumann graph is of a similar order of magnitude as predicted 
by our theory but its massive contribution is underestimated. More work needs 
to be done to capture the massive contributions for higher moments correctly. 
Our theory may be improved considerably by starting from a different exact ex- 
pression or by going beyond the Gaussian approximation in the massive modes. 
In the orthogonal case one should start from an expression that incorporates the 
symmetries of the wave function in all orders. 

3. The crossover regime, characterized 1 < a < 2 in ( 1260b . As limT?/ + the 
Gaussian Random Wave Model does not hold. In this regime we conjecture that 
the criteria ( 1256b or ( 1258b decide whether a sequence of graphs is asymptotically 
quantum ergodic or not. This conjecture for the crossover regime should be taken 
with much more care than the previous conjecture for the large mass regime. It 
does work for Neumann star graphs which have exponent a - \ and for which 
different methods revealed that asymptotic quantum ergodicity does not hold Q. 
These graphs have indeed a large number of masses with m ~ 1 jB. This number 
is of order 0{B) such that the limit in ( 1258b gives a constant. However, it has 
also been brought to our attention pTi that analogous criteria derived in Il22l l23l 
for the validity of Gaussian Random Matrix predictions for spectral correlation 
functions may lead to wrong conclusions for some borderline cases for which 
the analogous massive contributions are overestimated using the saddle-point ap- 
proximation to the corresponding exact variant of the supersymmetric cr-model. 
To some extent the prediction for the massive correction may be improved as 
outlined in the large mass regime. 

4. The non-universal small mass regime, a > 2 in ( 1260b . We conjecture that neither 
the Gaussian Random Wave Model nor asymptotic quantum ergodicity hold. In 
this regime the saddle-point analysis to the exact supersymmetric cr-model may 
break down completely. While it may not be trivial to prove this part of our 
conjecture rigorously our results give very strong evidence in favour of the con- 
jecture. 

8 Discussions 

Our main results are the formula ( 1250b for the autocorrelation functions C\a\ defined 
in (l42l l. and the formula ( 1255b for the fluctuations of an observable defined in 



55 



These formulae depend on the quantum graph only through the matrix R, and this 
matrix, which is defined in ( l35T l, only involves the underlying classical dynamics M. 
Hence, our results relate the statistical properties of the quantum energy eigenfunctions 
to properties of the classical dynamics on the graph. Moreover, they also reveal that 
the system dependency has no chance to vanish, and hence, a finite graph cannot be 
entirely described by the Gaussian Random Wave Model developed in Section [3?2l or 
even be quantum ergodic. These properties can only be met asymptotically in increas- 
ing sequences of graphs, that is in sequences of graphs where the number of bonds 
tends to infinity. In Section |7] classical criteria for such a sequence to be asymptoti- 
cally described by the Gaussian Random Wave Model or to be asymptotically quantum 
ergodic are formulated (the section concludes with a summary of the criteria and con- 
nected conjectures). The condition 7? — > for full universality, that is for the Gaussian 
model to be satisfied in the large graph limit, is more restrictive than the criteria for 
asymptotic quantum ergodicity. This is understandable since this latter property only 
depends on the second moment of the intensities, and the fluctuations 'Fv in ( |255t , 
which measure the deviation to ergodicity, can also decay in a non-universal way. 

The general formulae ( 12501 ) for the autocorrelation functions and, in particular, 
(I255l l for the fluctuations, have been obtained by a saddle-point analysis of the ex- 
act field-theoretical expression ( |148t . A comparison with the two periodic orbits ap- 
proaches in the subsections |3.5| and |4.2| reveals how the field-theoretical scheme exactly 
proceeds. The first term of 'Fv in ( |255t . which is universal, originates from our exact 
calculation on the zero mode manifold, and it coincides with the result predicted by the 
long diagonal orbits in Subsection [33] This draws a parallel between the zero mode, 
that is the uniform component of the classical map M, and long diagonal orbits. This is 
in fact not surprising since the zero mode is precisely the one that does not decay, and 
can thus survive in long orbits. The second term of Tv involves the system-dependent 
matrix R, that is the non-zero masses, and it coincides with the system-dependent con- 
tribution of the diagonal approximation exposed in Subsection 14.21 Hence, one de- 
duces that our field-theoretical approach discriminates between the diff'erent modes of 
the classical map M. The uniform component of M is treated in an exact way, which 
the diagonal approximation in 14.21 cannot do, while the massive decaying modes are 
treated in a perturbative way. 

It is also interesting to compare our results with those obtained by S. Gnutzmann 
and A. Altland in ll22l and 1231 concerning the asymptotic spectral two-point correla- 
tion function ^2(^) in a sequence of increasing quantum graphs. Their theory relates 
the function Riis) to the sequence of spectral gaps A of the matrices 1 - M. If the 
spectral gaps stay away from zero, the random matrix two-point correlation function is 
obtained in the limit of large graphs. Our condition /? — > for full universality requires 
all the non-zero eigenvalues of 1 - M to tend to one, which is obviously much stronger. 
Hence, even in situations where the Gaussian Random Wave Model does not hold, there 
is a possibility for random matrix theory to describe R2{s), but if the Gaussian Wave 
Model does hold, then R2is) must be universal. Moreover, if the sequence of spectral 
gaps vanishes as Am ~ B" as the number of bonds B becomes large, Gnutzmann and 
Altland's theory predicts different outcomes for R2{s) depending on the value of the 
positive number a. If o- < \,a random matrix behavior is reached, whereas a non-zero 
system-dependent contribution always remains if a > 1 . In the intermediate regime 
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a e [5, 1), the asymptotic two-point function R2(s) depends on the proportion of van- 
ishing modes, as explained in 1231 . In Subsection 17.21 we found that a < I imphes 
asymptotic quantum ergodicity, whereas a >2 generally forbids ergodicity. Therefore, 
universality for R2{s) implies asymptotic quantum ergodicity. However, in the domain 
a € [7,1), quantum ergodicity is always reached, whereas R2{s) can be non-universal. 
To conclude, let us mention some possible improvements of our method and some 
interesting directions for further research. In the main formula (I25Q| |, the system- 
dependent terms correspond to a Gaussian approximation around Q = cr^'^ in the direc- 
tions that are transverse to the saddle-point manifold. A true Gaussian approximation 
should expand the exact action to second order around every point of the saddle-point 
manifold. The correspondence between these two procedures has only been verified 
on the submanifold Z = Z^ with vanishing sources. Moreover and more importantly, 
in this expansion around the zero mode manifold, the higher order terms have not been 
controlled. Estimating these terms remains a major problem of this field-theoretical 
method. Note that only in the fully universal case /? — > one knows a small parame- 
ter {R itself) that one may use to order a systematic expansion. Besides, we have also 
shown that different but equivalent conventions in ( l76b lead to different outcomes by 
our second order expansion scheme. This implies that our second order expansion is 
not a systematic expansion in any intrinsic parameter of the quantum graph. Another 
question is whether the formula ( 12501 ) is suitable to describe other quantum systems 
if the matrix M is replaced with the Perron-Frobenius operator of a chaotic Hamilto- 
nian system. The field-theoretical method used here is probably difficult to generalize 
to other systems. An idea would be to develop a single periodic orbit approach that 
reproduces (1250b and transfer it to other types of systems. 

Acknowledgments This work has been partially supported by ORS AS and the Swiss 
Society of Friends of the Weizmann Institute of Science. The authors are grateful to 
Yan Fyodorov, Martin Sieber and Uzy Smilansky for their helpful comments. 

References 

[1] A. Backer, R. Schubert, R Stifter, Rhys. Rev. E 57, 5425 (1998). 

[2] A. Backer, R. Schubert, J. Phys. A 35, 527 (2002). 

[3] A. Backer, R. Schubert, J. Phys. A 35, 539 (2002). 

[4] F Barra, P Gaspard, J. Stat. Phys. 101, Nos. 1/2 (2000). 

[5] F A. Berezin, Introduction to superanalysis, D. Reidel Publishing Company, Dor- 
drecht, Holland, 1987. 

[6] G. Berkolaiko, Proc. Symp. Pure Math. 77, 247 (2009). 

[7] G. Berkolaiko, E.B. Bogomolny, J.R Keating, J. Phys. A 34, 335 (2001). 

[8] G. Berkolaiko, J.P Keating, J. Phys. A 32, 7827 (1999). 



57 



[9] G. Berkolaiko, J.P. Keating, U. Smilansky, Commun. Math. Phys. 273, 137 
(2007). 



[10 
[11 
[12 
[13 
[14 
[15 
[16 
[17 
[18 

[19 

[20; 
[21 
[22 
[23 
[24 
[25 
[26 

[27: 

[28 
[29; 
[30 
[31 
[32 
[33 
[34 



G. Berkolaiko, J.P. Keating, B. Winn, Phys. Rev. Lett. 91, 134103 (2003). 

G. Berkolaiko, J.P Keating, B. Winn, Commun. Math. Phys. 250, 259, (2004). 



G. Berkolaiko, B. Winn. larXiv:080 1.4 1041 /2, 2008. 

M.V. Berry, J. Phys. A 10, 2083 (1977). 

M.V. Berry, K.E. Mount, Rep. Prog. Phys. 35, 315 (1972). 

J. Bolte and S. Endres, Proc. Symp. Pure Math. 77, 247 (2009). 

A. Bouzouina, S. De Bievre, Commun. Math. Phys. 178, 83 (1996). 
Y. CoUn de Verdiere, Commun. Math. Phys. 102, 497 (1985). 

B. Eckhardt, S. Fishman, J.P. Keating, O. Agam, J. Main, K. Miiller, Phys. Rev. 
E 52, 5893 (1995). 

K. Efetov, Supersymmetry in disorder and chaos, Cambridge University Press, 
United Kingdom, 1997. 

Y.V. Fyodorov, A.D. MirHn, IJMPB 8, 3795 (1994). 

M. Feingold and A. Peres, Phys. Rev. A 34, 591 (1986). 

S. Gnutzmann, A. Atland, Phys. Rev. Lett. 93, 194101 (2004). 

S. Gnutzmann, A. Atland, Phys. Rev. E 72, 056215 (2005). 

S. Gnutzmann, J.P Keating, F Piotet, Phys. Rev. Lett 101, 264102 (2008). 

S. Gnutzmann, U. Smilansky, Adv. in Phys. 55, 527 (2006). 

S. Gnutzmann, U. Smilansky, J. Weber, Waves in Random Media 14, S61 (2004). 

M.C. Gutzwiller, J. Math. Phys 8, 1979 (1967). 

R.A. Horn, C.R. Johnson, Matrix Analysis, Cambridge University Press, 1990. 

S. Hortikar, M. Srednicki, Phys. Rev. Lett. 80, 1646 (1998). 

J.R Keating, Proc. Symp. Pure Math. 77, 279 (2009). 

J.P Keating, J. Marklof, B. Winn, Commun. Math. Phys. 241, 421, (2003). 

J.P Keating, S.D. Prado, Proc. R. Soc. Lond. A 457, 1855 (2001). 

V. Kostrykin, R. Schrader, J. Phys. A 32, 595 (1999) . 

T. Kottos, U. Smilansky, Phys. Rev. Lett. 79, 4794 (1997). 



58 



[35] T. Kottos, U. Smilansky, Annals of Physics 274, 76 (1999). 

[36] A.D. Mirlin, Phys. Rep. 326, 259 (2000). 

[37] F. Piotet, Statistical properties of the eigenfunctions on quantum graphs, PhD 
thesis, Bristol, 2009. 

[38] H. Schanz, T. Kottos, Phys. Rev. Lett. 90, 234101 (2003). 

[39] R. Schubert, Ann. Henri Poincare 7, 1085 (2006). 

[40] A.I. Shnirelman, Usp. Mat. Nauk 30, 265 (1975). 

[41] I. Smolyarenko, personal communication. 

[42] J.D. Urbina, K. Richter, J. Phys. A 36, 495 (2003). 

[43] J.D. Urbina, K. Richter, Phys. Rev. Lett. 97, 214101 (2006). 

[44] J.D. Urbina, K. Richter, Eur. Phys. J. 145, 255 (2007). 

[45] P. Walters, Ergodic theory: Introductory lectures. Lecture Notes in Mathematics 
458, Springer- Verlag, BerUn, 1975. 

[46] E.P Wigner, Group theory and its application to quantum mechanics of atomic 
spectra. Academic Press Inc, New York, 1959. 

[47] M. Wilkinson, J. Phys. A 21, 1173 (1988). 

[48] S. Zelditch, Commun. Math. Phys. 160, 81 (1994). 

[49] S. Zelditch, Commun. in Partial Differential Equations 19, Issue 9 & 10, 1565 
(1994). 

[50] S. Zelditch, M. Zworski, Commun. Math. Phys. 175, 673 (1996). 

[51] M.R. Zirnbauer, J. Phys. A 29, 7113 (1996). 



59 



